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Abstract. Motivated by models of fracture mechanics, this paper is devoted to the analysis of unilateral gradient flows 
of the Ambrosio-Tortorclli functional, where unilaterality comes from an irreversibility constraint on the fracture density. 
In the spirit of gradient flows in metric spaces, such evolutions arc defined in terms of curves of maximal unilateral slope, 
and are constructed by means of implicit Euler schemes. An asymptotic analysis in the Mumford-Shah regime is also 
carried out. It shows the convergence towards a generalized heat equation outside a time increasing crack set. 
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1. Introduction 

Many free discontinuity problems are variational in nature and involve two unknowns, a function u and a 
discontinuity set T across which u may jump. The most famous example is certainely the minimization of the 
Mumford-Shah (MS) functional introduced in [55] to approach image segmentation. It is defined by 



where VL C is a bounded open set, J^^ -1 is the (N — l)-dimensional Hausdorff measure, j3 > is a fidelity 
(constant) factor, and g 6 L°°(Q) stands for the grey level of the original image. In the resulting minimization 
process, we end up with a segmented image u : il \ T —> K and a set of contours r C fl. To efficiently tackle this 
problem, a weak formulation in the space of Special functions of Bounded Variation has been suggested and 
solved in [21], where the set T is replaced by the jump set J u of u. The new energy is defined for u £ SBV 2 (fl) 



where Vm is now intended to be the measure theoretic gradient of u. 

A related model based on Mumford-Shah type functionals has been introduced by Francfort & Marigo 
in [2j5] (see also [5]) to describe quasi-static crack propagation inside elastic bodies. It is a variational model 
relying on three fundamental principles: 

• the fractured body must stay in elastic equilibrium at each time (quasi-static hypothesis); 

• the crack can only grow (irreversibility constraint); 

• an energy balance holds. 

In the anti-plane setting, the equilibrium and irreversibility principles lead us to look for constrained critical 
points (or local minimizcrs) at each time of the Mumford-Shah functional, where u stands now for the scalar 
displacement while T is the crack. Unfortunately, there is no canonical notion of local minimality since the 
family of all admissible cracks is not endowed with a natural topology. The research of local minimizers of 
such energies has consequently become a great challenge, and a lot of works in that direction have considered 
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global minimizers instead, see [17114124] . In the discrete setting, one looks at each time step for a pair (ui,Ti) 
minimizing 

(«,r)^l/ \vu\ 2 dx + je N -\r), 

1 Jn\r 

among all cracks T D Ti-i and all displacements u : il \ T — > R satisfying an updated boundary condition, 
where Ti—i is the crack found at the previous time step. A first attempt to local minimization has been carried 
out in [18] where a variant of this model is considered. At each time step the L 2 (f2)-distance to the previous 
displacement is penalized. More precisely, denoting by Ui-i the displacement at the previous time step, one 
looks for minimizers of 

(u,r)->i/ \Vu\ 2 dx + .yf N - 1 {Y) + \\\u-u^ 1 \\l Hn)l (1.2) 
1 Jn\r 

on the same class of competitors than before, where A > is a fixed parameter. We emphasize that this 
formulation only involves some kind of local minimality with respect to the displacement. A notion of stability 
which implies another local minimality criterion has been introduced in [32| . It focuses on what the author calls 
"accessibility between two states" . In the case of global minimization, when passing from one discrete time to 
the next, all states are accessible. From the point of view of [18], a state u is accessible from if and only if 
there is a certain gradient flow beginning at which approaches u in the long-time limit. The main idea in 
[32] is that a state u is accessible from if and only if both states can be connected through a continuous 
path for which the total energy is never increased more than a fixed amount. 

While static free discontinuity problems start to be well understood, many questions remain open concerning 
their evolutionary version. Apart from the quasi-static case, the closest evolution problem to statics consists in 
finding a steepest gradient descent of the energy, and thus in solving a gradient flow type equation. A major 
difficulty in this setting is to define a suitable notion of gradient since the functional is neither regular nor 
convex, and standard theories such as maximal monotone operators |11] do not apply. To overcome the use of 
differential or subdiffcrential, a general theory of gradient flows in metric spaces has been introduced in [35] 
where the notion of gradient is replaced by the notion of slope, and the standard gradient flow equation is recast 
in terms of curves of maximal slope (see [1] for a detailed description of thissubject). The construction of such 
generalized evolution is usually performed by means of an implicit Euler scheme (corresponding to a discrete- 
in-time steepest descent) whose limits are refered to as De GiORGi's minimizing movements, see [1120] , The 
minimizing movements of the Mumford-Shah functional have been first considered in [2], and further developed 
in [12] . Motivated by the crack growth model, the authors actually apply the iterative scheme with respect to the 
variable u while minimizing the energy with respect to T under the constraint of irreversibility. More precisely, 
they minimize energy (|1.2[) exactly as in [18] with A replaced by (2<5) _1 , where S > is the time step. Showing 
compactness of the resulting discrete evolution as 5 — > 0, they obtain existence of "unilateral" minimizing 
movements of the Mumford-Shah energy (we add here the adjective unilateral to underline the irreversibility 
constraint on the evolution). In any space dimension, the limiting displacement u(t) satisfies some kind of heat 
equation (in a very weak sense), and an energy inequality with respect to the initial time holds. Assuming that 
admissibe cracks are compact and connected, they improve the result in two dimensions showing that u(t) solves 
a true heat equation in a fractured space-time domain, and that the energy inequality holds between arbitrary 
times. With the (probable) aim of relating the unilateral minimizing movements of |2|12] to some generalized 
gradient flow, a notion of unilateral slope for the Mumford-Shah functional has been introduced in [19]. Up to 
our knowledge, no precise relation between these objects has been found yet. 

The Mumford-Shah functional enjoys good variational approximation properties by means of regular energies. 
Constructing gradient flows for these regularized energies and taking the limit in the approximation parameter 
could be a way to derive a generalized gradient flow for MS. It was actually the path followed in [30] where 
a gradient flow equation for the one-dimensional Mumford-Shah functional is obtained as a limit of ordinary 
differential equations derived from a non-local approximation of MS. Many other approximations arc available, 
and the most famous one is certainely the Ambrosio-Tortorelli functional defined for (u,p) € [i? 1 (fi)] 2 by 

AT e (u, p) := 1 J^{n e + p 2 )\Vu\ 2 dx + \J Q ( £ I V H 2 + \( l ~ P?) dx + f jf (« ~ 9? dx . 

The idea is to replace the discontinuity set T by a (diffuse) phase field variable, denoted by p : fl — !• [0,1], which 
is "smooth" and identically in a e-ncighborhood of T. Such energies are of great importance for numerical 
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simulations in imaging or brittle fracture, see |7l8j . From the mechanical point of view, it is interpreted as a 
non-local damage approximation of fracture models, where p represents a fracture density. The approximation 
result of |5l6j (see also [55]) states that AT e T-converges as e — > to MS (in the form with respect to a 

suitable topology. For the static problem, it implies the convergence of AT e -minimizers towards MS-minimizers 
by standard results from r-convergence theory. However, the convergence of general critical points is a priori 
not guaranteed. Positive results in this direction have been obtained in |25l34j for the one-dimensional case. 
Note that it is also of numerical interest to investigate the behavior of general critical points of AT e . Indeed, 
while existence of minimizers is easy to prove, it is not straightforward to compute them numerically since 
the functional AT £ is not globally convex (but only separately convex). The algorithm used in [5J consists in 
performing alternate minimizations in each variable, and it only provides critical points of AT £ . 

The Ambrosio-Tortorelli approximation of quasi-static crack evolution is considered in [29], where the ir- 
reversibility constraint translates into the decrease monotonicity of the phase field t i— y p(t). The main result 
of [29] concerns the convergence of this regularized model towards the original one in [24]. Motivated by the 
formulation of a model of fracture dynamics, a hyperbolic evolution related to the Ambrosio-Tortorelli func- 
tional is also studied in [33], but the asymptotic behavior of solutions as e — > is left open. A first step in that 
direction is made in [15] where the analysis of a wave equation on a domain with growing cracks is performed. 
Concerning parabolic type evolutions, a full gradient flow of the Ambrosio-Tortorelli functional is numerically 
investigated in [27| for the purposes of image segmentation and inpainting. 

The object of the present article is to present a new notion of unilateral gradient flow for the Ambrosio- 
Tortorelli functional taking into account the irreversibility constraint on the phase field variable. Motivated 
by [19] , we adopt the general framework of gradient flows in metric spaces [3] , defined here as curves of maximal 
unilateral slope (see Definition 13. 7[) . Existence is obtained through the minimizing movements method, where 
the discrete Euler scheme is precisely the Ambrosio-Tortorelli regularization of the one studied in |2ll2j . As 
in [53], the irreversibility of the process is encoded into the decrease monotonicity of the phase field variable, 
and leads to constrained minimization problems. At the discrete time level, given an initial data (uo,po), one 
recursively defines pairs (v,i, pi) by minimizing at each time U ~ iS, 

(u,p) t-^AT e (u,p) + —\\u-u i - 1 \\l a ^ n)) 0-3) 

among all u and p < p%—i, where (ui—i,pi—i) is a pair found at the previous time step. The objective is 
then to pass to the limit as the time step 6 tends to 0. A main difficulty is to deal with the asymptotics of 
the obstacle problems in the p variable. It is known that such problems are not stable with respect to weak 
i/ 1 (r2)-convergence, and that "strange terms" of capacitary type may appear |13ll6j . However, having uniform 
convergence of obstacles would be enough to rule out this situation. For that reason, instead of AT ei we consider a 
modified Ambrosio-Tortorelli functional with p-growth in Vp with p > N. By the Sobolev Imbedding Theorem, 
with such a functional in hand, uniform convergence on the p variable is now ensured. We define for every 
(u,p) E H 1 ^) x iy 1 'P(fi), 



S s {u,p):=- / (r k + p i )\Vu\'dx + / \v P \r + —\l-p\'>\dx + Z / {u-gfdx, 

2 Jn Jn\ V pe J 2 J n 

where a > is a suitable normalizing factor defined in (|2.2[) . Note that an immediate adaptation of [35] shows 
that £ £ is still an approximation of MS in the sense of r-convergence. 

Considering the incremental scheme (| 1 . 3[) with £ e instead of AT E , we prove that the discrete evolutions 
strongly converge as 6 — > to continous evolutions t t— > (u £ (t), p e (t)) that we call unilateral minimizing move- 
ments (see Definition 13. 4[) . These limiting evolutions turn out to be I/ 2 (0)-steepest descents of £ e with respect 
to u in the direction of non- increasing p's, i.e., curves of maximal unilateral slope in the spirit of [4119] (see 
Theorems 15.11 & 15. 2\ . The Euler-Lagrangc equation for u £ is given by 



d t u e - div((?7 e + p 2 e )Vu e ) + P(u E - g) = in Vt x (0, +oo), 

OH 

— - = on dtt x (0,+oo), 

av 



(1.4) 



while the irreversibility and minimality conditions for p e are 
1 1— > p £ (t) is non-increasing , 

(1.5) 

£ e (u e (t), p e (t)) < £ e (u e (t),p) for every t > and p <E W ltP (tt) such that p < p e (t) in fl. 
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The system (|1.4[) - (|1.5p is supplemented with the initial condition 

(u £ (0),p £ (0)) = (u o ,p ) inO. 
In addition, we prove that the bulk and diffuse surface energies, defined by 

\ J (Ve + p e {t) 2 )\Vu £ {t)\ 2 dx 

and 

t ^ f (—\ Vpe{t )\r + ^-\l- Pe (t)Adx 
Jn V P P £ J 

are respectively non-increasing and non-decreasing, a fact which is meaningful from the mechanical point of 

view. Moreover, the total energy is non-increasing, and it satisfies the following inequality: for a.e. s £ [0, +00) 

and every t > s, 

e e (u s (t),p e (t))+ f \\d t u e (r)\\l HQ) dr<S e (u s (s),p £ (s)). 

J s 

Note that the inequality above is reminiscent of gradient flow type equations, and that it usually reduces 
to equality whenever the flow is regular enough. In any case, an energy equality would be equivalent to the 
absolute continuity in time of the total energy. The reverse inequality might be obtained through an abstract 
infinite-dimensional chain-rule formula in the spirit of |37j . In our case, if we formally differentiate in time the 
total energy, we obtain 

j t £ £ (u e (t),p £ (t)) = (d u £ £ (u £ (t),p £ (t)),d t u £ (t)) + (d p £ £ (u £ (t),p e (t)),d t p £ (t)) . 
From (|1.5[) we could expect that 

(d p £ £ (u £ (t), p £ (t)), d t p £ {t)) = , (1.6) 

which would lead, together with \1A\ , to the energy equality. Now observe that (| 1 .6[) is precisely the regu- 
larized version of Griffith's criterion stating that a crack evolves if and only if the release of bulk energy is 
compensated by the increase of surface energy (see e.g. Section 2.1]). Unfortunately, the chain-rule above 
is not available since we do not have enough control on the time regularity of p £ . In the quasi-static case, one 
observes discontinuous time evolutions for the surface energy. Since the evolution law for p £ is quite similar to 
the quasi-static case (see [29]), we also expect here time discontinuities for the diffuse surface energy. Adding 
a parabolic regularization in p could be a way to improve the time regularity, but unfortunately, it would also 
break the increase monotonicity of the surface energy. 

A natural continuation to the qualitative analysis of Ambrosio-Tortorelli minimizing movements is now to 
understand their limiting behavior as e — > 0. We stress that the general theory on T-convergence of gradient 
flows as presented in [38 315] does not apply here since it requires a well defined gradient structure for the T- 
limit. A specific analysis thus seems to be necessary. In doing so, we show that (u £ ,p £ ) tends to (u*, 1) for some 
mapping 1 1— > u*(t) taking values in S BV 2 (O) , and solving (in a weak sense) the equation 

{<9 t u* — div(Vw*) + f3(u* — g) = in £1 x (0, +00) , 
Vu*-v = on dfl x (0, +00) , (1.7) 

u*(0) = u . 

In addition, using the monotonicity of the diffuse surface energy, we are able to pass to the limit in (|1.5|) . It 
yields the existence of a non-decreasing family of rectifiable subsets {T(t)} t >o of Q such that J u ,(t) C T(i) for 
every t > 0, and for which the following energy inequality holds at any time: 

\ I \Vu t (t)\ 2 dx + Jf N - 1 (T(t)) + ^ f (u4t)~g) 2 dx+ f \K( S )\\ 2 LHn) ds 

<o / |Vu | 2 cfe + ^ / {u -g) 2 dx. 
1 Jn 1 Jn 

Comparing our result with [12] , we find that solves the same generalized heat equation with an improvement 
in the energy inequality where an increasing family of cracks appears. The optimality of this inequality and the 
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convergence of energies remain open problems. Note that the convergence of the bulk energy usually follows by 
taking the solution as test function in the equation. In our case it asks the question wether SBV 2 (£l) functions 
whose jump set is contained in J Ut ( t ) can be used in the variational formulation of (| 1 . T[) . It would yield a weak 
form of the relation 

((u+(t)-„-(t))*^ = onr(t), 

where uf(t) denote the one-sided traces of u*(t) on T(t). This is indeed the missing equation to complement (|1.7[) . 
and it is intimately related to the finiteness of the unilateral slope of the Mumford-Shah functional evaluated 
at (ii*(t),r(t)) (see [TH Proposition 1.1], and Section HT3l) . 

The paper is organized as follows. In Section^ wc define in details the Ambrosio-Tortorelli and Mumford- 
Shah functionals as well as the functional setting of the problem. In the first part of Section [31 we introduce 
two implicit Eulcr schemes generating two different types of unilateral minimizing movements. The first one is 
based on global minimization, while the second one relies on alternate minimization as in |7l8j . In the second 
part of Section [3J we define and analyse the notions of unilateral slope and curve of maximal unilateral slope 
for the Ambrosio-Tortorelli functional. In Section [4l we establish compactness of discrete evolutions when 
the time step tends to zero. In order to get sharp energy inequalities, we are led to extend the concept of 
De Giorgi's interpolation to our unilateral setting. Then we study in Section [5] some qualitative properties 
of unilateral minimizing movements, showing in particular that they are curves of maximal unilateral slope. 
Finally, Section [6] is devoted to the asymptotic analysis as e — > 0, where we show that unilateral minimizing 
movements of the Ambrosio-Tortorelli functional converge to solutions of the generalized heat equation (|1.7I) 
outside time-increasing cracks. 

Notations. For an open set U C M. N , we denote by ^(U; R m ) the space of all finite R m -valucd Radon measures 
on U, i.e., the topological dual of the space %)(U;M. m ) of all R m -valucd continuous functions vanishing on dU. 
For m = 1 we simply write ^(U). The Lebesgue measure in M. N is denoted by , while Jtf' 1 ^^ 1 stands for 
the (N — l)-dimcnsional Hausdorff measure. If B\ is the open unit ball in R*, we write ujn '■= J? N (Bi). We 
use the notations C and = for inclusions or equalities between sets up to J^^" ^negligible sets. For two real 
numbers a and 6, we denote by a A 6 and a V b the minimum and maximum value between a and b, respectively, 
and a + := a V 0. 



2. Preliminaries 

Absolutely continuous functions. Throughout the paper, we consider the integration theory for Banach 
space valued functions in the sense of Bochner. All standard definitions and results we shall use can be found 
in [11] Appendix] (see also [23]). We just recall here some basic facts. If X denotes a Banach space, we say 
that a mapping u : [0,+oo) — > X is absolutely continuous, and we write u G AC ([0, +oo); X) , if there exists 
m G L 1 (0, +oo) such that 

\\u(s) — u(t) \\x < J m(r) dr for every t > s > . (2-1) 

If the space X turns out to be reflexive, then any map u G AC([0, +oo); X) is (strongly) derivable almost 
everywhere. More precisely, for a.e. t G (0, +oo), there exists u'(t) G X such that 

u(t) — u(s) ... . 

> u (t) strongly m X as s — > t . 

Moreover v! G L (0, +oo; X), vl coincides with distributional derivative of u, and the Fundamental Theorem of 
Calculus holds, i.e., 



u(t) — u(s) 



J u (r) dr for every t > s > . 



If further the function m in (|2.1[) belongs to L 2 (0,+oo), then we write u G AC 2 ([0, +oo); X), and in that case 
we have v! G L 2 (0, +00; X). 

Special functions of bounded variation. For an open set U C M. N , we denote by BV(U) the space of 
functions of bounded variation, i.e., the space of all functions u G L l (U) whose distributional gradient Du 
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belongs to ^(U;M. N ). We shall also consider the subspace SBV(U) of special functions of bounded variation 
made of functions u G BV(U) whose derivative Du can be decomposed as 

Du = Vu3? N + o+ - u-)v u ^e N - x L J u . 

In the previous expression, V« is the Radon-Nikodym derivative of Du with respect to Ji? N , and it is called 
approximate gradient of u. The Borel set J u is the (approximate) jump set of u. It is a countably Jt? 1 ^ -1 - 
rectifiable subset of U oriented by the (normal) direction of jump i/ u : J u — > § w_1 , and u ± are the one-sided 
approximate limits of u on J u according to v u , see [3]. We say that a measurable set E has finite perimeter in U 
if xe G BV(U), and we denote by d*E its reduced boundary. We also denote by GSBV(U) the space of all 
measurable functions u : U —> M such that (— M V u) A M £ SBV(U) for all M > 0. Again, we refer to [3] for 
an exhaustive treatment on the subject. Finally we define the spaces 

SBV 2 (U) := {u G 5W(C/)flL 2 ([/) : Vu G L 2 (U;R N ) and J^" 1 ^) < oo} , 

and 

GSBV 2 (U) := {u G GSBV{U) n L 2 (C/) : Vu G L 2 ({/;M JV ) and 3^ N ~ X {J U ) < oo} . 

Note that, according to the chain rule formula for real valued BV-functions, we have the inclusion SBV 2 (U) n 
L°°(U) c GSBV 2 (U) (see e.g. Theorem 3.99]). 

The following proposition will be very useful to derive a lower estimate for the Ambrosio-Tortorelli functional. 
It is a direct consequence of the proof of [51 Theorem 10.6] (see [TUl Theorem 16] for the original proof). 

Proposition 2.1. Let C be a bounded open set, let {u„}„ 6 n C H 1 ^!) fl L°°(Cl) be such that 

sup ngN ||tt n ||ioo(o) < o°; an d let {E n } ne tq be a sequence of subsets of £1 of finite perimeter in f2 such that 
l Jrf? N ~ 1 (d*E n fl fl) < oo. Assume that u n — > u strongly in L 2 (Jl), and that J£ N (E n ) — > 0. Set- 



SU Pn(EN ■ 

ting u r , 



(1 — XE n )u n G SBV (fi) fl L°°(p,), assume in addition that sup ngN ||Vu 



n\\L 2 (n;R N ) 



< oo. Then 



u G SBV 2 (n)nL°°{Q), and 



u n — > u strongly in L 2 (Q) , 

u n u weakly* in L°°(Q) , 

Vm„ Vu weakly in L 2 (Q;R N ) , 

2Jt? N - 1 (J u ) < HmMJf N - 1 (d*E n n fi) . 



The Ambrosio-Tortorelli & Mumford-Shah junctionals. Let us now consider a bounded open subset fl 
of M. N , p > N, f3 > 0, and g G L°°(f2). For e > and G (0, 1), we define the Ambrosio-Tortorelli functional 
£ E : L 2 (fl) x U>(fl) ->• [0, +oo] by 



1 



(% + p 2 )\Vu\ 2 dx + 



£ —\Vp\r + ^\l- P r 
n \ P PE 



£ £ {u,p) := < 



+ 2 / (u - gf dx 



if (u,p) G 1 (£1) x W x *{n) , 



otherwise , 

where p' := p/(p — 1) and a is the normalizing factor given by 

»=-(' I* 1 -j** *r -(if- 

The Mumford-Shah functional £ : L 2 (S1) — > [0, +oo] is in turn defined by 

I [ \\7u\ 2 dx + .J*t? N ~ 1 (J u ) + ^ [ (u-g) 2 dx if u G GSBV 2 (fl) , 

otherwise . 



£{u) :-- 



-OC 



(2.2) 



(2.3) 



It is well known by now that the Ambrosio-Tortorelli functional approximates as e — > the Mumford-Shah 
functional in the sense of T- convergence, as stated in the following result, see |5l28j . Let us mention that 
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Theorem 12.21 is not precisely a direct consequence of |5l28j . In [5], the case p = 2 is adressed, while [55] deals 
with energies having the same p-growth in Vu and Vp (recall that p > N > 2). However, a careful inspection 
of the proof of [5SJ Theorem 3.1] shows that the T-convergence result still holds for £ £ . 

Theorem 2.2. Assume that r/ £ = o(e). Then £ e Y -converges as e — >• (with respect to the strong L 2 (to) x L p (to)- 
topology) to the functional £q defined by 



£ (u,p) 



£{u) ifue GSBV 2 (to) and p = l into, 
+oo otherwise . 



3. Unilateral minimizing movements & curves of maximal unilateral slope 
3.1. The discrete evolution schemes 

Throughout the paper, we shall say that a sequence of time steps d := {<5 4 }i 6 N» is a partition of [0, +oo) if 

S* > , sup Si < +oo , and > <T = +oo . 
»>i jri 

To a partition 5 we associate the sequence of discrete times {i 4 }i<=N given by t° := 0, t % := J2]=i & f° r * — 1) 
and we define the time step length by 

|5| := sup 5 l . 

»>i 

To an initial datum uo £ H 1 ^) n L°°(to), we shall always associate (for simplicity) the initial state p% 
determined by 

Pq ■= argmin £ e (u ,p). (3.1) 

It is standard to check that the above minimization problem has a unique solution (by coercivity and strict 
convexity of the functional £ £ (uq, •)), and it follows by minimality that < p\ < 1. Given a partition S of 
[0, +oo), we now introduce two discrete evolution schemes starting from (uq,pI). 

Scheme 1 (global minimization): Set (u°,p°) := (uq,Pq), and select recursively for all integer i > 1, 

(«*,/>*) G argminj^Kp) + ± \\u - u 1 ' 1 1|| 2(0) : («, p) G H x (to) x M /1 ' p (fi) , p < p l ~ x m flj . (3.2) 

Scheme 2 (alternate minimization): Set (u°,p°) := (uq,/Oo), and define recursively for all integer i > 1, 



gmin {^(u.p*- 1 ) + i ||«- ^-^^(n) : « £ H 1 ^)} , 
{f e « p):pe W^ p (ft) , p < p 1 - 1 m to} . 



p := argmm 

While it is straightforward to check that the minimization problems in Scheme 2 admit (unique) solutions, 
the well-poscdncss of Scheme 1 requires some little care. Since the sublevel sets of £ e are clearly relatively 
compact for the sequential weak 77 1 (fi) x W 1,p (Q)-topology, one may apply the Direct Method of Calculus 
of Variations to solve (|3.2|) . We only need to show that the constraint in (|3.2| is closed, and that £ £ is lower 
semicontinuous with respect to weak convergence. 

Lemma 3.1. Let {(«„, p n )}neN C .f/ 1 (SI) x W 1,p (to) be such that (u n , p n ) — (u, p) weakly in H 1 (to) xW 1 ^ (to). 
Then, 

£e(u,p) < hminf £ s (u n ,p n ) ■ (3.3) 

n— >oo 

Moreover, if for each n G N, p n < p in to for some p G W 1,p (to), then p < p in to. Finally, assuming that 
£ e (u n , p n ) — > £ s (u, p) as n — > oo, then (u n , p n ) — > (u,p) strongly in H 1 ^) x W 1,p (to). 
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Proof. Step 1. The sequence {(u n , p n )} being weakly convergent, it is bounded in x W 1,p (fl). Therefore 

p n — > p in ^ (fi) by the Sobolev Imbedding Theorem. Hence p < p in fl whenever p n < p in ft for every n £ N. 
Then p n Vu n pVit weakly in L 2 (Q), and consequently, 



/ (r] £ + p 2 )\Vu\ 2 dx < lim inf / (r) E + pl)\Vu n \ 2 dx . 
Jn n ^>°° Jn 



Since all other terms in £ e are clearly lower semicontinuous with respect to the weak convergence in x 
W^ifl), we have proved Q . 



Step 2. Let us now assume that £ £ (u n ,p„) — > £ e {u 1 p). We first claim that 

(% + p 2 )\Vu\ 2 dx = lim / (r? £ + p 2 n )\Vu n \ 2 dx . (3.4) 
i n ^°°Jn 

Indeed, assume by contradiction that for a subsequence {rij} we have 

/ (?7e + p 2 )\Vu\ 2 dx < lim inf / + p 2 n . )\ Vu n | 2 dx . 

Jn -7^°° Jn 3 

Using the fact that u n — > u strongly in L 2 (tt), we deduce from Step 1 that 

lim SeiUn.jPm) 
3->oo 

>liminfi / (r lE + p 2 l .)\Vu n] \ 2 dx + \immi [ ( S —\Vp \P + -2L|1 - Prlj A dx + f / (u - g) 2 dx 
j^oo 2 J n J j^oo J n \ p p'e 'J 2 J n 

> £e(u,p) , 

which is impossible. Therefore (|3.4[) holds. Then, combining the convergence of £ e (u n , p n ) with (|3.4[) . we deduce 
that ||pn||w 1 'P(n) —> \\p\\w 1 -p(n), whence the strong lU^^f^-convergence of p n . 

It now remains to show that u n —¥ u strongly in H l (Vl). Using the uniform convergence of p n established in 
Step 1, we first estimate 

P 2 - Pn\\^u n \ 2 dx < (su-p\\\7u k \\ L 2 ia . M N ) )\\p 2 - pl\\ Lx{n) — > 0. 



Then we infer from (|3.4|) that 

{rie + p 2 )\Vu n \ 2 dx = f + p 2 n )\Vu n \ 2 dx + I (p 2 ~ p 2 n )\\7u n \ 2 dx — »■ / (r/ £ + p 2 )|Vu| 2 cfe . 

Consequently ||wn||_f/ 1 (n) —> \\u\\h 1 (Q), whence the strong ff 1 (0)-convergence of u n . □ 

Remark 3.2. Let us now briefly comment both algorithms. The first scheme is very much in the spirit of classical 
minimizing movements, and it is the analogue of the algorithm introduced in [12] for studying minimizing 
movements of the Mumford-Shah functional constrained to irreversible crack growth. Scheme 1 can be seen as a 
regularization of the one in [12] since the Ambrosio-Tortorelli functional L-converges as e — > to the Mumford- 
Shah functional (see 028]). However, a practical drawback of Scheme 1 is that the pair (u l ,p l ) obtained at 
each time step might not be unique since £ e is not strictly convex (although is is separately strictly convex). 
This lack of uniqueness may generate some troubles from the point of view of numerical approximations. For 
that reason, it is of importance to consider an alternate scheme which uniquely defines the pair {u 1 , p 1 ). It turns 
out that Scheme 2 is actually the algorithm used in numerical experiments for quasi-static evolution in brittle 
fracture (see |7|8j ). We shall prove that both schemes give rise to the same time continuous model, once the 
time step length \8\ is sent to zero. 

We state below some important pointwisc estimates on the iterates {(u*, p l )}i 6 N which easily follow from 
minimality and standard truncation arguments. 

Lemma 3.3. Let {(«\ /0 I )}«eN be a sequence obtained from either Scheme 1 or Scheme 2. Then, for every i £ N, 
l|tt*IU-(n) < max{||« |U-(n), IMU~(fi)} and < p l+1 < p i < 1 in 0. (3.5) 
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From a sequence of iterates {(u l ,p I )}jgN obtained from cither Scheme 1 or Scheme 2 and a partition S 
of [0, +oo), we associate a discrete trajectory (us,ps) : [0,+oo) — > x M /1,p (f2) defined as the Ze/t 

continuous piecewise constant interpolation of the (w l ,p z )'s below. More precisely, we set 

us(Q) = u Q , p s (0) = p E , 

and for t > 0, 

M ):=u * ifteC^.f]. (3.6) 

By analogy with the standard notion of minimizing movements, we now introduce the following definition. 

Definition 3.4 (Unilateral - Alternate - Minimizing Movements). Let uq £ i7 1 (17)nL°°(r2). We say that 
a pair (u, p) : [0, +oo) — > L 2 (Q) x L p (f2) is a (generalized) unilateral minimizing movement (resp. a (generalized) 
unilateral alternate minimizing movement) for £ e starting from (uo,/9q) if there exist a sequence {£fc}fcgN of 
partitions of [0, +oo) satisfying \Sk\ ^ 0, and associated discrete trajectories {(us k , ps k )}kefi obtained from 
Scheme 1 (resp. Scheme 2) such that 

(u Sk (t),p Sk (t)) — > (u(t),p(t)) strongly in L 2 {fl) x L p (fi) for every t > 0. 

k— ► oo 

We denote by GUMM(uq, p§) (resp. GUAMM(uo, pg)) the collection of all (generalized) unilateral minimizing 
movements (resp. (generalized) unilateral alternate minimizing movements) for £ e starting from (uq,Pq). 

Remark 3.5. At this stage we do not claim that the collections GUMM(uq, p§) and GUAMM(uq, Pq) are 
not empty. This will be proved in the next section through a compactness result on discrete trajectories (see 
Lemmas rOlfe 14. Ill and Corollary [4. 14j) . 



3.2. Curves of maximal unilateral slope 

In the spirit of [3], we introduce in this section the notion of L 2 (f2)-unilatcral gradient flow for the Ambrosio- 
Tortorelli functional in terms of curves of maximal unilateral slope, accounting for the quasi-stationnarity and 
the decrease monotonicity constraint on the phase field variable p. To this aim, we first define the unilateral 
slope of £ e which is analogous to the one introduced in [TO] for the Mumford-Shah functional. 

Definition 3.6. The unilateral slope of £ e at (u, p) € i? x (0) x W 1,p (f2) is defined by 

\ac u ^ v I {£e{u,p) - £ £ (v,p)) + ! , 1 

\o£ s \(u,p) := hmsup sup < n n — : p G W ,P (S1), p< pm.il} , 

v^u p { \\v -u||i,2(n) J 

where the convergence v — > u holds in L 2 (Jl). The functional \d£ e \ is then extended to L 2 (Jl) x L p (il) by setting 
\d£ s \(u,p) := +oo for (u,p) £ ff 1 ^) x W 1 '^). 

The unilateral slope being defined, we can now define curves of maximal unilateral slope for the Ambrosio- 
Tortorelli functional. 

Definition 3.7. We say that a pair (u,p) : (a, 6) — > L 2 (f2) x L p (0) is a curve of maximal unilateral slope for 
£ e if u 6 AC 2 (a, 6; L 2 (f2)), p is non-increasing, and if there exists a non-increasing function A : (a, b) — > [0, +oo) 
such that for a.e. t £ (a, 6), £ e (u(t), p(t)) = X(t), and 

A'(t) < -\\\u\t)\\ 2 L2[n) - \\d£ e \ 2 (u{t),p{t)) . (3.7) 

Our definition of curve of maximal unilateral slope is motivated by the following proposition which parallels 
[1 Theorem 1.2.5]. 

Proposition 3.8. If (u, p) : (a, b) —> L 2 (Q) x L p (fl) is a curve of maximal unilateral slope for £ e , then 

\\u'(t)\\ LHQ) = \d£ £ \(u(t),p(t)) for a.e. t E (a, b) . (3.8) 
Moreover, if t i— > £ £ (u(t), p(t)) is absolutely continuous on (a,b), then 

£ e (u(t), p{t)) + / \\u'(r)\\ L 2fQ\dr = £ £ (u(s),p(s)^ for every s and t <G (a, 6) with s < t . 

J s 
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Proof. Let A be as in Definition 13.71 Since A is non- increasing, A has finite pointwise variation in (a, b). Let us 
consider the set 

A := [t 6 (a, b) : £ B (u(t), p(t)) = \(t), A and u are derivable at i} , 

and observe that Jz? 1 ^, b) \ A) = 0. 

Let t <E A. Since A is non- increasing, we have X(t) < 0, and thus 

|A'(*)| = -X'(t) = Inn m^m = Um £Mt),pit))-£s{u{ 8 ), P (s)) 
sit, sGA s — t sit, seA s — t 

Using the fact that p(s) < p(t) when s > t (by the non-increasing property of t i— > p(t)) and the strong 
L 2 (f2)-continuity of u, we infer that 

IVmi <- r (£ £ (u(t),p(t)) - £ £ (u(s),p)) + \\u(s) - u(t)\\ LHn) 

A i < hnisup sup ^ r—— — '- < \d£ £ \{u(t), p(t)) u (i) L 2 (u) ■ 

slts£Ap<p(t) \\U(S) - U(t)\\ L 2 {n) S-t 

On the other hand, \\'{t)\ > \\\u'{t)\\ 2 L2[n) + \\d£ £ \ 2 (u(t), p{t)) by (02]), and fllTS]) follows as well as the fact 
that \'{t) = -\\u'{t)\\l, m . 

Finally, if i i — >• £ £ (u(t), p(t)) is absolutely continuous on (a, 6), then for every s,t € (a, b) with s < t, 

£ £ (u(t), p(t)) - £ £ (u(s),p(s)) = j A'(r) dr = - j \\u'(r)\\ 2 L2{n} dr , 
which completes the proof of the proposition. □ 

We state below necessary and sufficient conditions for the finiteness of the slope, as well as an explicit formula 
to represent it. This is one of the milestone of our entire analysis. 

Proposition 3.9. Assume that Q, has a 'tf 1,1 -boundary , and let D[\d£ £ \) be the proper domain of \d£ e \. Then, 

D{\d£ e \) = \ {u,p) e H 2 {tt) x W U} (VL) : ^ =0 in i? 1/2 (<9ft), and 
\ av 

£e(u, p) < £ £ (u, p) for all p e W 1 ^^) such that p < p in Q j . (3.9) 

In addition, for (it, p) G D(\d£ £ \), 

\d£ £ \(u,p) = ||div((Tfe + P 2 )Vu) - 0(u - g)\\ L2m , (3.10) 

and 

\M\n»((i) < C £ (l + \\V p\\ L p(n-M")) a (\d£ £ \(u, p) + /3\\u -g\\ L 2 (n} + \\u\\ H i {n) ) , 
where a 6 N is the smallest integer larger than or equal to pj (p — N), and C £ only depends on r\ £ , p, N , and £1. 
The proof of Proposition ^. 91 is based on the following auxiliary regularity result. 

Lemma 3.10. Assume that Q has a e ta 1 ' 1 -boundary. For f <G L 2 (0) and p € W 1,p (fl), let u 6 i/ 1 (fi) be a 
solution of 

' -div((r? e + p 2 )Vu) = f in H- l {Vl) 

(3.11) 

(?7e + p 2 )Vu ■ v = m J ff- 1 / 2 (9fi). 

du 

Then u E H 2 (n), — = in H l / 2 (dVL), and 
av 

< C e (l + ||V / 9|| iP(a . K ^ ) ) a (||/|| i 2 (n) + IMIhmo)) , 
where a GN is as in Provosition \3.9\, and C £ only depends on r\ £ , p, N, and f2. 



Proof. Step 1. We claim that 

du 



-0 in H-^ 2 {dfl). (3.12) 
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To prove this claim, we first rewrite the equation as 

rVp ■ Vu H 

Ve + P 2 Ve + P 



Au = - 7 2 Vp ■ Vu + - , ^ 2 in^'(fi). (3.13) 



Hence Am G L q (il) with g := 2p/(p+ 2) by Holder's inequality. Then we observe that q' := q/(q — 1) < 2* since 
p > TV, so that (51) by the Sobolev Imbedding. Hence the linear mapping 



ipeH 1 ^)^ / (Vu • Vip + (Au)ip) dx 
Jn 



is well defined and continuous. Consequently, u admits a (weak) normal derivative on dft which belongs to 
the dual space H~ 1 / 2 (dCl), and for any ip £ i/ 1 (fl), 



du 

at/ / (H- l /2(ao),Hi/2(an)) Jn 



/ (Vu • Vip + (Au)(p) dx 
Jn 

[ (ri e + p 2 )Vu ■ V [ — ^—5-) dx+ [ ( Au H ^rVp- Vu | (pete. 

in \Ve + P 2 J Jn\ Ve + P 2 J 



We observe that in the second equality above, we have used the fact that G H 1 ^) whenever 1/3 G i? 1 (il). 

Indeed, 



V A _ V<p 2p<pVp ^ r2 



\Ve+ P J Ve + P (Ve + P ) 

since ip G i 2 *(fi), p G L°°(fl), and Vp G L p (f2) with p > N. In view of ([3~TTj) we have 

f (r h + p 2 )\7u-v(—?—) dx= f ~^dx, 
Jn \Ve + P 2 J JnVe + P 2 

and by (|STTS|> . 

^—^Vp ■ Vu ) ipdx = — [ — — da; , 

% + P / % + P~ 

from which (|3.12[) follows. 

Step 2. We now prove that u G H 2 (Q). By the previous step, u G i? 1 (fi) satisfies 

Au G L«(fi) , 
|^ = 0in J ff- 1 /2(^). 

By elliptic regularity (see e.g. [31] Proposition 2.5.2.3 & Theorem 2.3.3.6]), we deduce that u G W 2 ' q °(fl) with 
<lo -=q= ^f2, and 

IMIm/ 2 -™^) < C(l|Au|| £TO (n) + ||w|U«o(n)) . 

for some constant C > only depending on N, p, and f2. Observing that the function t h-> t/(r] e +t 2 ) is bounded, 
we derive from (|3.13[) and Holder's inequality that 

IMIw 2 '<Jo(n) < C e (||/||i2(n) + ||Vp|| L p(o ; k« ) ||Vu|| L 2 ( j 2 . R ]v ) + ||u|| L 2 (f2) ) 

< C e (l + \\Wp\\ LH a. RN) )(\\f\\ L 2 iu) + ||u|| H i ( n)) , 

where we used the fact that qo < 2. By the Sobolev Imbedding, we have u G W 1,q ° (CI), and thus Vu- Vp G L qi (fi) 
with 

1 1 1 2Np 
■ , i-e., q% ' 



qi P q* (N-2)p + AN 

Note that q\ > 2 if and only if p > 2N, so we have to distinguish the case p > 27V from the case p < 2N. 
Case 1). Let us first assume that p > 2N. Then Vu • Vp G L 2 (Cl) with 

|jVu • Vp|| L 2 (n) < C||Vu • Vp|| Lgi(0) < C||Vp|| L p (n . R iv ) ||Vu|j i ,* (0;KJV) < C||V/o||jrp( n . I ,/f)||u|| H «>, eo (n) ■ 
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Using again (|3.12[) - (|3.13[) and the elliptic regularity, we infer that u £ H 2 (fl) with the estimate 

IMIff 2 (f2) < C(\\Au\\ L 2 {n) + \\u\\ L 2 {n) ) 

< Ce{\\f\\L*(n) + ||Vp|| L p( O ;R«)||u|| W /2,, 0(o) + ||u||i,2(n)) 

< C e {\ + ||V P || LP(W ) 2 (||/|| i2(n) + ||«|| ff i ( n)) ■ 

Case 2). If p < 2N then qi < 2, and we have u £ W 2 ' qi (fl) by (|3TT^ - (|3"Tl"3"|> and elliptic regularity, with the 
estimate 



\\u\\w 2 -n(<n) < C E (\\f\\ L 2 {n) + \\Vp\\Lp(n)\\u\\w 2 'm(n) + ||«IU 2 (n)) 

< C e (l + \\V P \\ LPm ) 2 (\\f\\ L 2 m + \\u\\ H i {a) ) . (3.14) 

In particular, Vw, £ L qi (fi) by the Sobolev Imbedding since q\ < 2 < N. We then continue the process by 
setting 

1 1 1 2Np 
i.e., qi ■ 



qi ' p " *" (JV-2i)p + 2(i+l)JV 

as long as q^i < 2, that is i < a. Since q a -i > 2, iterating estimates of the form Q3.14p we obtain 

\\u\\H 2 (n) < C e (\\f\\ L 2 (n} + ||Vp|| LP(J2;R iv ) ||u|| H/ 2,, c ,_ 2(n) + \\u\\ L 2 {n) ) 

< c s (i + ||Vp|U, (W ) a (||/|U 2(fi) + \\u\\ m(u) ) , 

and the proof is complete. □ 

Proof of Proposition [3791 Step 1. Let us consider a pair (u,p) such that \d£ e \(u,p) < oo. For ip £ i? 1 (il) 
with ip ^ 0, we estimate 

\d£ e \(u,p)>lun S np £Au ' p) ~f {u ~ S(p ' p) > ^ ( [ (r? £ + p 2 )Vu . V*> + /3 / (u-^ds). (3.15) 

By density of 7J 1 (fi) in L 2 (0) and the Riesz representation Theorem in L 2 (il), we deduce that there exists 
/ £ L 2 (Vt) such that 

(r] E + p 2 ) Vu ■ Vip dx + (3 I (u — g)ipdx= I fipdx 
Jfi Jn 

for all <p £ H 1 ^). Hence u solves (|3.11[) with / = f — j3{u — g). We then infer from Lemma f3 . 101 that u £ H 2 (il), 
and that ^ = in H 1 ^ 2 (dfl). Next, taking ip £ i? 1 (fi) such that ||y||x, 2 (S7) = 1j integrating by parts in (|3.15[) . 
and passing to the supremum over all such tp's yields the lower bound 

\d£ s \(u,p) > \\div(( % + p 2 )Vu) - f3(u - g)\\ L2{n) . 

We now claim that the following minimality property for p holds: 

£ e (u, p) < £ £ (u, p) for all p £ W l - p (VL) such that p < p in O . (3.16) 

Since \d£ e \(u 1 p) < +oo we can find sequences {v n } C i/ 1 (0) and {p n } C W 1,P (Q) such that v n — > u strongly 
in L 2 (n), 

p n = argmin {£ e (t>„, p) : p £ W 1,p (f2), p < p in 17} for each n £ N , 

£ s (Vn, Pn) < £e{u, p) , (3-17) 

and 

£ e (u,p) - £ s (v n ,p n ) 
hmsup 1, n < \d£ e \[u,p). (3.18) 

n-s-oo \\V n — U\\ L 2(q) 

By (|3.17p the sequence {Vw„} is uniformly bounded in L 2 (£l;M. N ). Hence, for a suitable subsequence (not 
relabeled), 

\Vv n \ 2 ^f N \_n -± \Vu\ 2 ^ N \_ n + p 
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weakly* in J?(R N ) for some nonncgative Radon measure p G JV(R N ) supported in tt. Let us now consider the 
following functionals on W 1,P (SY) defined by 

_ ... (£ e (v n ,p) Hp<p, . \£ e (u,p) + \ [ (Ve + P 2 )dp ifp<p, 

■r n (p) := < and F(p) := < ^ Jn 

[+oo otherwise, |^ +00 otherwise. 

Note that by the Sobolev Imbedding W 1,p (f2) ^ ^°(f2), the functional T is well defined on the space W 1,p (£l). 

We claim that T n T-converges to T for the sequential weak W 1,p (f2)-topology. Indeed, the upper bound 
is immediate since T n (p) — > J-(p) for each p G W l ' p {Vt). For what concerns the lower bound, if {p„} C 
W 1,P {£1) is such that liminf„ .F ra (p„) < oo and p„ — p weakly in W 1,p (£l), then for a subsequence {nk} we 
have limfc T nk (p nfc ) = lim inf n T n (p n ) , and p nk — > p in ^°(fi) by the compact imbedding iy 1,p (f2) =-> ^°(f2). 
Consequently p < p in fi, and 



/ #jV«nJ a ds-> / / 5 2 |Vu| 2 d.T+ /_p 2 dp. 



Since the remaining terms in the energy T n arc independent of n and lower scmicontinuous for the weak 
W 1,p (f2)-convergence, we deduce that 

F{p) < lim inf J"„(p„) , 

and the T-convergence is proved. 

Since the sublevel sets of J- n are relatively compact for the sequential weak T4 /1,p (r2)-topology (uniformly 
in n), we infer from the T-convergence of T n towards JF that 

£ £ (v n ,p n ) = min T n -> min T. 

On the other hand, by (|3.17|) and (|3.18|) we have £ £ (v n ,p n ) — > £ e (u, p) from which we deduce that 
£ £ (u, p) = min T = min {Tip) : p G VF 1,p (f2), p < p in fi} . 

We conclude from this last relation that p, = and that (|3.16[) holds. 

Step 2. Conversely, we show that if a pair (u, p) belongs to the set in the right hand side of (|3.9[) , then 
\d£ E \(u,p) < oo and \d£ £ \(u,p) < ||div((?7 £ + p 2 )Vu) - (3(u - g)\\ L2{ny 

Consider a pair (it, p) G ff 2 (ft) x W^fi) satisfying |^ = in H x l 2 (dn) and 

£ £ (u,p) < £ £ (u,p) 

for all p G W /1 ' p (^) such that p < p in D,. Note that u G W 1 ^^) for every r < 2* by the Sobolev Imbedding, 
and since p > N, the product V« • Vp belongs to L 2 (f2) and p G L°°(J7). Hence, 

div((?7 £ + p 2 )Vu) = (r) £ + p 2 )Ait + 2pVp • Vu G L 2 (0), 

and consequently, it is enough to check that 

\d£ e \(u,p) < ||div((r7 £ +p 2 )V U )-/3( U - 5 )|| i2(0) . 

Consider a sequence {v n } C i? 1 (fi) converging strongly to w in L 2 (f2) such that 

, 1f |, v ,. / {£e{u,p) -£ e {v n ,p)) + . . . n \ 

\o£ £ \(u, p) = lim sup < - - - ' — : p G W ' P (U), p < p m it > , 



\v n - u\\L 2 (n) 



and let 



Then 



so that 



p n = argmin {£ £ (v n , p) : p G W 1,p (fi) such that p < p in fi} 



sup (£ e (u,p) - £ £ (v n ,p)) < (£ £ (u,p) - £ £ {v n ,p n )) 
p<p 



\»C \l ^ V (£e{u, p) - £ £ {v n , p n )) 

\d£ £ \(u,p) = lim ^ '—. (3.19) 
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If for infinitely many ro's we have £ s (v n , p„) > £ £ (u, p), then |9£ e |(it, p) = and there is nothing to prove. Hence 
we can assume without loss of generality that £ £ (v n ,p n ) < £ £ (u,p). In particular, {p n } is uniformly bounded 
in W 1,P (Q), and {v n } is uniformly bounded in H 1 ^). As a consequence, for a subsequence v n — u weakly in 
H 1 ^) and p n p* weakly in W 1,p (f2). From Lemma |3. II we infer that p* < p in f2, and 

£ £ (u,p*) < liminf £ e (v n ,pn) < limsup£ E (v n , p n ) < £ s (u,p) . (3.20) 

n ^oo n— >oo 

By the minimality property of p, we have that £ £ (u, p) < £ £ (u 1 p„) which leads to £ £ (u, p) = £ £ (u,p*). By 
uniqueness of the minimizcr (due to the strict convexity of £ £ (u 1 ■)), we deduce that p* = p. Then Lemma |3. II 
and (|3.20p with p* = p shows that p n — »• p strongly in W 1 ' P (Q). 
We now estimate 

£ £ {u,p) - £ £ (v n ,p n ) < £ £ (u,p n ) - £ s (v n ,p n ) 

< / {"He + Pn)^ u ■ (Vu — Vw„) dx + P I (u — g)(u — v n ) dx 
Jn Jn 

= - / (u - u„)(div((r? £ + p n )Vu) - /3{u - g)) dx . 



Note that in the last equality, there is no boundary term since = in H 1 ^ 2 (dfl). Moreover, since u € H 2 (Q) 
and p n € W 1,p (fi), we have div((ry e + p 2 )Vii) € L 2 (f2). Applying Cauchy-Schwarz Inequality we obtain 

£,.(-», p) - £ £ {v n ,p n ) 2 

ii [I < div (t? £ + p n )Vu) - f3(u - g)\\ L 2(n) . 

\\v n ~ u|| L 2 (a) 

Since H 2 (Vt) ^ W 1,r (fi) for every r < 2* and p n -> p strongly in M /1,p (0), we get that 

div((? ?£ + p 2 n )Vu) = (r) E + p 2 )Au + 2p n Vp n ■ Vu — > (r/ e + p 2 )Au + 2pVp • Vu = div((?7 £ + p 2 )Vu) 

n— >oo 

strongly in £ 2 (f£). Hence 

lim £ ^>P)~ £ ^n,Pn) < ||diy(( + p2)Vu) _ (u _ }|| 
n->oo \\v n - u|| L 2 (0) 

Together with (|3.19|) . this last estimate gives the desired upper bound for the slope \d£ £ \(u, p). 

Step 3. Let (u, p) G D(\d£ £ \). By the previous steps, / := — div((?/ £ + p 2 )Vu) + [3(u — g) e £ 2 (S1) and u solves 
(|3.1ip with / = / — (3(u — g). Applying Lemma [3.101 we hnd that 

h\\ H *(n) < C £ {1 + ||Vp|| L p (f2;R « ) ) Q (||/|| i 2 (0) + \\u\\ H i (n) ) 

< C £ (l + \\Vp\\ LP{ n-M-)) a (\d£ s \(u, p) +P\\u- g\\ L 2 {Q) + \\u\\ m{n) ) , 
and the proof is complete. □ 



The expression of the slope and the characterization of its domain provided by Proposition [379] enables one 
to show the lower semicontinuity of \d£ £ \ along sequences with uniformly bounded energy 

Proposition 3.11. Assume that £1 has a 1,1 -boundary. Let {(u n , p n )}neft C L 2 (fl) x L p (fl) be such that 
sup ngN £ £ {u n , p n ) < oo and {u ni p n ) { u -, P) strongly in L 2 (fl) x L p (fl). Then, 

| d£ £ | (u, p) < lim inf | d£ £ \ (u n , p n ) ■ 

n— >oo 

Proof. Let us assume without loss of generality that lim inf „ \d£ £ \(u n , p n ) < oo, and extract a subsequence 
{n k } such that 

liminf \d£ £ \(u ni p n ) = lim \d£ £ \(u nkl p nk ) ■ 

n— >oo k— s-oo 

Since £ £ {u nk , p nk ) is uniformly bounded with respect to k , we deduce that the sequence { ( Wn fc , p nk ) } is uniformly 
bounded in i? 1 (fi) x W 1,p {Vt). Moreover (u nk ,p nk ) € D(\d£ £ \), and as a consequence of Proposition 13. 9\ we 
deduce that {u nk } is uniformly bounded in H 2 (n), and that dl t ) ^ k = in H 1 / 2 (dVL). Whence p Uk p weakly 
in W 1,P (Q), u nk — ^ u weakly in H 2 (Q) for a (not relabeled) subsequence, and |^ = in H 1 / 2 (dfl). By the 
Sobolcv Imbedding we get that p nk — s- p in ^f°(f2), while u nk — > u strongly in 1 (SI) . Thanks to the uniform 
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convergence of p nk to p, we may argue as in the proof of Proposition 13.91 Step 1 with minor modifications, to 
show that the functionals Tk ■ W 1,p (il) —> [0, +oo] defined by 

F k (p) := < . (3.21) 

I +oo otherwise , 

F-converge (with respect to the sequential weak W 1,p (fi)-topology) to the functional F : W 1,p (fi) — > [0, +oo] 
given by 

Hp)--={ ee{U,P) if/5 " P ' (3-22) 
+oo otherwise . 



Indeed, the lower bound inequality can be obtained as in the proof of Proposition 13.91 Step 1, while the upper 
bound requires the following argument. Given p £ W 1,P (Q) satisfying p < p in ft, we set ps := p — S where S > 
is small. Since p nk — > p uniformly in f2, we have ps < Pn k in ^ whenever k large enough (depending only on S). 
Hence, 

lim limsup F k (ps) < F(p) , 

s l° fc->oo 

and we obtain from {ps}s>a & suitable recovery sequence for p through a diagonalization argument. 
Since 

p Uk = argmin F k (p) , 

and Pn k p weakly in W ,p (Cl), we infer from the T-convergence of F k toward F that 

p = argmin J- (ft) ■ 

By the expression of the domain of the slope provided by Proposition ^. 91 we infer that (u,p) <G D(\d£ E \). From 
the established convergences of (u nk ,p nk ) we deduce that 

div((?7 e + p 2 nk )Vu nk ) = (% + pl k )Au nk + 2p nk Vp nk ■ Vu„ fc 

(f) e + p 2 )Au + 2pVp ■ Vu = div((r/ £ + p 2 )Vu) 
weakly in L 2 (Q). Using now the expression of the slope given by Proposition 13.91 we conclude 

liminf \d£ e \(u n ,p n ) = lim \d£ e \{u nk , p nk ) 

n—too k— >oo 

= lim ||div((? ?e + pi JVu„J - (3(u nk - g)\\ L 2 {n) 

k—>oo 

> \\div(( % +p 2 )X7u)-f3(u-g)\\ L 2 m = \d£ e \(u,p), 
which ends the proof. □ 

We next prove that the energy is continuous along convergent sequences with uniformly bounded slope. 
Proposition 3.12. Assume that has a % 71,1 -boundary. Let {(u n , p n )}n£fi C L 2 (f2) x L p (Vl) be such that 

sup{5 e (M n ,/9„) + \d£ s \(u n , p n )} < oo, 

and (u ni p n ) — ► (u, p) strongly in L 2 (fl) x L p (fl). Then £ £ (u n ,p n ) — > £ £ (u,p) as n — > oo. 

Proof. Arguing as in the proof of Proposition 13.111 we have u n — 1 u weakly in H 2 (fl) and p n — 1 p weakly in 
W 1,p (n) with (u,p) e D(\d£ e \). By the Sobolev Imbedding, p n -> p in and u n -> u strongly in ff^O). 

Hence the functional F n : W hp (Q) -> [0, +oo] defined by (|3T2T]) (with n in place of rife) T-converges (with 
respect to the sequential weak W 1,p (fi)-topology) to the functional F : W lp (il) — > [0, +oo] given by (|3.22[) . By 
the convergence of the minimum values, we infer that 

£ e (u n ,p n )= min F n — > min F = £ £ (u, p) , 
and the proposition is proved. □ 
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4. Compactness of discrete trajectories 

The purpose of this section is to obtain compactness properties of discrete- in-time evolutions {(ug k , ps k )}keN 
associated to a sequence {5fc}fcgN of partitions of [0, +00) as the time step length \8 k \ tends to zero. In a first 
subsection we shall define a suitable notion of Dc Giorgi interpolants in order to derive from Scheme 1 a precise 
discrete energy identity (see (|4.1ip ). This interpolation procedure extends the standard De Giorgi interpolation 
as defined in [3] to our unilateral setting. For Scheme 2, a direct computation shall provide an optimal energy 
inequality (see (|4.12jl ). From these energy estimates and from the regularity result obtained in Proposition 13.91 
(concerning the slope of the Ambrosio-Tortorelli functional), we will derive a priori estimates leading to the 
compactness of {(us k , ps k )}keN in various functional spaces. 

Throughout this section we fix an arbitrary uq G H 1 (S})PiL°° (fi) , and we consider the function determined 
by (|3.1[) . A sequence {<5fc}fc £ N of partitions of [0, +00) satisfying \8k\ ~ > is also given. We write 

i 

8 k =: {4} ieN . , t% := , and t< := £ 5{ for i > 1 . 

3=1 

For each fceNwe consider a discrete trajectory (uk, Pk) ■ [0, +00) — > H 1 ^) x W 1,p (il) associated to 8k which 
is obtained from either Scheme 1 or Scheme 2 (and (|3.6I) ). We also assume that every element of the resulting 
sequence {(life, p/c)}fceN arises from the same scheme as in Definition 13.41 To simplify the notation, we write 

«, P i):= M4),^(4))- 

We next define for every k G N a further piecewise constant interpolation p k : [0, +00) —> W lp (tt) from the 
iterates {p},}^ setting p k (0) = pg, and for t > 0, 

PkW^pr 1 i{t€(ti-\ti). (4.1) 

We also consider the piecewise affine interpolation v k ■ [0, +00) — > _ff 1 (f2) of the u^'s defined for each k G N by 

v k {t) := u*" 1 + t —^« - if * 6 t**" 1 .*!] ■ ( 4 - 2 ) 

4.1. De Giorgi interpolants for Scheme 1 

We present in this subsection a notion of unilateral De Giorgi interpolation very much in the spirit of [4j 
Chapter 3]. Up to minor modifications due to the unilateral constraint, the proofs of the results below follow 
closely their analogues in [H Chapter 3]. For clarity reasons we have decided to present all the details. 

Let us fix k G N, and assume that the iterates {(u l k , p k )}ieN are obtained from Scheme 1. For an integer 
i > 1 and t > 0, we define ^, fe (r; • ) : if 1 (Q) x W^ity -> [0, +00] by 

^ k (r;u,p) := S e (u,p) + ^rllw-wr^lla^) • 
Then we consider the function <pi >k : (0, +00) — > [0, +00) given by 

0i lfc (r) :=inf {$ a (r;u, p) : (u,p) G H^Q) x W^(fi) , p < p^ 1 in fl} . 

To each (pi,k(r) we associate the (resolvent) set ^im{t) c i/ 1 (SI) x W 1,p {^t) defined as 

J l .k(r) := argmin^ )fc (r; •) . 

We observe that J?i,k{j) ^ by the argument used in the proof of Lemma |3~T1 and that (u k , p k ) G fli,k{$k)- 

Definition 4.1 (De Giorgi interpolation). We consider (uk, p k ) '■ [0, +00) — > H l {VL) x W 1,p (il) an arbitrary 
interpolation of the iterates {(u k , p l k )}ien (obtained from Scheme 1) satisfying 

M4)>M4)) = (4,4) for*GN, 

(u k (t),Pk(t)) G Mr) ^t = ^ +re(f k \V k ). 

Let us emphasize that, at a formal level, a De Giorgi interpolation corresponds to a unilateral minimal energy 
path connecting the discrete values {{u k , p k )}i^. To investigate the main properties of (u k ,p k ) it therefore of 
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importance to look for the regularity of the energy along such interpolation. To this purpose it is useful to 
introduce the functions d^ fe : (0, +00) — > [0, +00) defined by 

d+ fc (r) := sup \\u T - u\~ l \\ L 2 [n) and dr fc (r) := inf \\u T - u^ 1 \\ L 2 (n) . 
(■u T , PT )ej?i, h (T) (.^,p r )eJi, k {T) 

First notice that 

\dt(T)\<C, (4.3) 



for some constant C > independent of k, i, and r. Indeed, we can argue as in Lemma 13731 to obtain the uniform 
bound 

IK||l«>(Q) < max{||M ||L-(o), I|5||l~(o)} for every r > and every (u T ,p T ) € c /i,k( T ) > 
from which estimate (14.31) follows. 



Lemma 4.2. Let i > 1, < n < T2, and (u Tj ,p Tj ) € a fi,k{'Tj) f or j = 1)2. VFe /iai/e 

ft <Mr 2 ) < ^(n) < ^k- 1 ,^- 1 ); 

(Mj ||u Tl -w^ 1 || i 2 (fi) < ||w T2 - w^ 1 || L 2 (n) ; 
fmj d+ fe (r 1 )<d- fc (r 2 )<d+ fc (r 2 ). 
In addition, there exists an (at most) countable set Mi. k C (0, +00) smc/i i/iai 

d tk( T ) = d i,k( T ) f° r ever V T e (°> +°°) \ M.fc • (4-4) 

Finally, 

hm0 iife (r) = 4K" 1 ^r 1 )- ( 4 - 5 ) 

Proof. It is straightforward to check (%), which shows that r 1— > (t) is non-increasing. Then, by minimality 
of (u T .,p T ,), we easily estimate 

« — 1 II 2 II., »,»-l||2 



^ + 5 e (« ri ,p T J < h£ £ (u T2 ,p T2 ) 



= T, \\ u r 2 -U k ||i2 (a) + 0i,fc(T2) < — \\U T2 -U k || i 2 (a) + +£ e (u Tl ,p Tl ) , 

and then derive 

T 2 - T l I, t-l||2 ^ T 2 - r l || t-l||2 

IL=(n)<-^ \\ u T 2 -u k || L2(n) , 

ZT1T2 ^TlT 2 

whence (ii) and (m). 

As a consequence of (jz), the functions r 1— > d~ fc (r) are non-decreasing, and therefore continuous outside an 
(at most) countable set Afi,k C (0, +00). Then, for r G (0, +00) \ A/i,jt, 

dr fc (r)<d^(r)<limdr fc ( S ) = dr fc (r), 

and (|4.4[) is proved. 

Let us now consider an arbitrary sequence r n 4- 0, and (u Tn ,p Tn ) £ flifrijn)- From ("i) we first deduce that 
£r(«r n ,Pr n ) < Ss^ 1 J^) and || Ut „ - ^H;^ < 2r n £ £ (u[7 1 , p^) -* . 

Hence we can find a subsequence (not relabeled) such that (u Tn , p Tj ,) K _1 ,p*) weakly in if 1 (ft) x W 1>p (n) 
for some 6 W 1,p (fi). From Lemma T3. II we derive that p* < p\~ x in J7, and 

^e(«fc -1 )P*) < liminf S £ (u Tn ,p Tn ) < Uminf iifc (r„) < limsup </> ijfe ( r «) < Zei^l' 1 , p^T 1 ) ■ 



n— foo n— »oo 



On the other hand, we have £ e (ul ,p*) > £e{ul iP\ *) by minimality f i^ anc j ^j^g (j) ik (j n ) — > 
£ E {u % k ~ ,p\~ ) as n — > cxd. Then (|4.5[) follows through the standard argument on the uniqueness of the limit. □ 

Remark 4.3. It follows from the previous lemma that the functions r n- df k (T) are Borel measurable. 
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dtjr) 



We next intend to obtain a discrete energy equality involving the term — acting as a "discrete time 
derivative" of u k . 



Proposition 4.4. For each integer i > 1, the function r i— > 4>i,k(j) is locally Lipschitz in (0, +oo), and 



»t,fc 



(r) 



« fe (-)) 2 (d^(r)) 5 



dr vw 2r 2 
in particular, for every {u T ,p T ) G ^^(t), 



2r 2 



£ e (u r ,p r ) + — ||tir-itfc 1 Hi2 (f2) + i 



and /or eac/i i > 1, 



i 



Proof. Let us fix arbitrary < t\ < t-2 and (u Tj ,p Tj ) € ^i,k{jj) f° r 3 = 1>2- We first observe that 



ds 



1 



for every t G (0, +oo) \7Vi,fe 



d ( 7 = f £ ( M r 1 ,pr 1 )- 







/o 


a 





d+» 


/ 

/o 


a 



(4.6) 



•7) 



(4.8) 



0t,k(n) - 4 ,fc(T 2 ) < $i,fc(ri;u T2 ,p T2 ) - $ iife (r 2 ;u T2 ,/9 T2 ) 



2tit 2 



■|l W r 2 - Ufc 1 |Il2(o) 



Similarly, we obtain that 



(f>i,k(Tl) - 0i,fc(T 2 ) > 



T2 - n 

2riT 2 



llL2(n) 



From the arbitrariness of (u Tj , p Tj ) we infer that 



(d+ fe (n)) 2 < ^ k ( Tl ) - ^ k (r 2 ) < (d" fc (r 2 )) 2 



2tiT2 T2 — Ti 2riT2 

which shows, together with (|4.3p . that r ^ k (r) is locally Lipschitz in (0, +oo). For r G (0, +oo) \ A/^fc, we 
can pass to the limit as t\ — > r and T2 — > r in the inequality above to derive (|4.6j) . 
Then, for < r < r, integrating (|4.6[) between ro and t yields 

2 



4>i,k{r) + 



1 



f 




'TO 


(7 



dcr = c/>i,k(ro) 



In view of (|4.5|) . it now remains to let tq 4- to recover (|4.7j) . Finally for each integer i > 1, we have (u l k , p k ) G 
^i,k{^i), and so identity (|4?7|> yields (|4~8|) . □ 



We end this subsection by estimating the slope of the Ambrosio-Tortorelli functional at the De Giorgi 
interpolants in terms of the "discrete time derivative" of u k . 

Lemma 4.5. For an integer i > 1, and t > 0, let (u Tl p T ) E J"i,k{j). Then (u T ,p T ) G DQd£ e \), and 

\d£ s \(u T ,p T ) < -\\u r - Mr 1 HL2(S2) ■ 
In particular, for every t > we have {u k (t),p k (t)) G D(\d£ £ \), and 

\d£ e \{u k {t),~Pk{t)) < G k (t) 
where G k : (0, +oo) — > [0, +oo) is defined by 

<k(r) 



(4.9) 



G k {t) := 



ift = tr 1 + r€(ti- \ty 



(4.10) 



Proof. Let (v,p) G x W 1,p (Cl) be such that p < p\ in fi. Then, by minimality of (u T ,p T ) 



£ e (u T ,Pr) ~£ e {v,p) < —\\v 

It 



l k IIl 2 ( o) 



4 1 !lz, 2 (fi) 



<2^(H»-»A 



L2(n) + IK - u k \\mn)) \\v - u T \\ L 2 {n) 
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Dividing this inequality by \\v — u r \\ L 2^ and taking the supremum over p, we derive that 

(£ e {u Tl p T ) -£ e (v,p)) + (£ £ {u T ,p T ) -£ e {v lP )) + 1 ._ x ._ x 

SU P n n < sup n 1, <-^-{\\v-u k \\l»(si) + \\ut - u k \\mn)) ■ 

P< P r h-u T \\ L 2 (n) p < p i-i \\v-u T \\ LHn) 2r 

Taking the limsup as v — > u T in L 2 (f2) then yields the desired result. □ 
4.2. A priori estimates and energy inequalities 

We first state a priori estimates based on discrete energy (in)equalities, which will later be essential to get 
compactness properties of the discrete trajectories. 

Lemma 4.6. If the discrete trajectory (u k , Pk) is obtained from Scheme 1, then for every i G N, 

1 /•** , . 1 



£e(ui,pi) + -J o \\v' k {s)\\l, m ds+-^ \Gk(s)\ 2 d s = £ E (u ,p £ ). (4.11) 

Proof. In view of the definition (|4. 10[) of Gk, iterating (|4.8I) from j = 1 to j = i yields (|4.11l) . □ 

For Scheme 2, we only obtain an (optimal) energy inequality between two arbitrary discrete times. 
Lemma 4.7. If the discrete trajectory (u k , Pk) is obtained from Scheme 2, then for every integers j > i, 

. . A 

^K.Pfc) + / Ktolli^O) ds ^ £ e( U lPk) ■ ( 4 - 12 ) 
Proof. By minimality, for each I € {i + 1, . . . , j}, u k satisfies 



^dx + I (r) e + (pi^fWui -V^dx + 13 j (u% - g)vdx = 



in d k Jn Jn 

for all ip € iJ 1 (fi). In view of this equation and the minimality of p k , we first compute 

£ e (u k -\ Pk - 1 ) = £ e (ulpir 1 )+ j ( % + (pir 1 f)Vui-V(u k ~ 1 -u e k )dx 

+ p I (4 - g^- 1 - u{) dx + \( ( % + (pir'n^K- 1 4)\ 2 <** + § / Kr 1 - 4) 2 ■ 

In L Jn z Jn 



u 



and then estimate 



2 rtt 



ui — u. 



dx>£ E (ui,pi)+ KOOIll^n)^- ( 4 - 13 ) 
J t e k 

Iterating this last inequality from I = i + 1 to I — j yields (|4.12l) . □ 



Proposition 13.91 then provides an additional H 2 (Q) estimate uniform in k for both schemes. 

Proposition 4.8. Assume that ft has a % n ' 1 -boundary. Then u k belongs to L 2 OC ([0, +oo); H 2 (fl)), = in 
L 2 (Q,+oo]H 1 /' 2 {d^l)), and for every T > 0, 

u k (t)\\ 2 HHn) dt<C E (T+l), (4.14) 



n 



for some constant C £ > independent of k. 

Proof. Step 1. We first consider the case where u k is obtain from Scheme 1. Let t > such that t s (i^ -1 ,^) 
for some integer i > 1. By minimality Uk{t) solves 

i-dW((n e +p 2 (t))\7uk(t)) = -v' k (t)-(3(uk(t)-g) intf- 1 ^), 
\(»/e + Pfc(t))V«i k (t)-i/ = mH-V*(dn). 

From Lemma l3~T0l we deduce that u k {t) £ H 2 (n), = in H x / 2 {dn), and the estimate 

IN(*)||ff*(n) < Ce(l + l|Vpfc(i)|| L P(j 2;R iv ) ) a (||ufe(i)|| L 2 (fi) +/3||ii fc (t) -5||z,2 (n) + \\u k {t)\\ H i (n) ) . 
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In view of (|4.11|) and Lemma 13.31 we deduce that 

\\u k (t)\\m(o) <C e (K(*)lli»(n) + l)> 
for some constant C £ > independent of fc, and (|4. 14[) follows integrating this last inequality between and T. 

Step 2. In the case where Uk is obtained from Scheme 2, we simply repeat the argument of Step 1 replacing pk 
by , and using (|4.12j) instead of (|4.1ip . □ 

4.3. Compactness 

We are now in position to to establish some compactness results for the sequences {w/cjfceN and {pfcjfceN- We 
start with {pfc}fc £ N and {p^}keN- 

Lemma 4.9. There exist a subsequence k n — > oo and a strongly measurable map p e : [0, +oo) — > W lp (fl) such 
that p kn {t) -± p e {t) weakly in W 1,p (fi) for every £ > 0. In addition, p e S L°°(0, +oo; W rl ' p (0)) ) p e (0) = pg, and 

< < Ps( s ) < 1 m 17 /or every t > s > . 

Proof. By Lemma f3.3[ pfc : [0, +oo) — > L 1 (fi) is monotone non- increasing, and < pk(t) < 1 in f2 for every 

1 > 0. By a generalized version of Holly's selection principle (see [35j Theorem 3.2]), we deduce that there exists 
a subsequence k n — ► oo and a map p e : [0, +oo) — > L 1 (fi) such that Pk n (t) — 1 p e (i) weakly in L 1 (fi) for every 
i > 0. On the other hand, since 

£s{u ,Po) < £ e {u 0l l) < |]Vu ]|L2 (n . K iV ) + ^\\u - g\\ 2 L 2 {n) , 
we derive from Lemmas 14.61 & 14.71 that 

sup \\pk(t)\\ w i, P{n ) < C £ , 
t>o 

for some constant C £ > independent of k. Therefore, Pk n (t) p £ (t) weakly in W 1,p (£l), and Pk n (t) — > Peit) 
in ^ (il) for every t > by the Sobolev Imbedding Theorem. In particular p e (t) £ W 1,p (f2) for every t > 0, 
and by lower semicontinuity, 

SUplIPe^)!!^ 1 -^^) < C e . 
t>0 

Moreover, since < p k {t) < Pfc(s) < 1 in f2 whenever s < t, we deduce from the uniform convergence that 
< Pe{t) < pe(s) < 1 in n for every t > s > 0. 

Since p £ : [0, +oo) — > W lp (ft) is a pointwise weak limit of a sequence of measurable (locally) simple 
functions, we deduce that p £ : [0, +oo) — >• W 1,p (Vt) is weakly measurable, hence strongly measurable thanks to 
the separability of W 1,p (tt) and Pettis Theorem. □ 

Through the same argument we obtain the convergence of the sequence {p/T}fceN (defined in (|4.1|) ). 

Lemma 4.10. Let {k n } n em be the subsequence given by Lemma \4-9\ There exist a further subsequence (not 
relabeled) and a strongly measurable map p~ : [0,+oo) — > W 1,P (Q.) such that p^~ (t) — p~(t) weakly in W ,p (Cl) 
for every t > 0. In addition, p~ g L°°{0, +oo; W lfP (n)), p~(0) = p%, and < p £ (t) < p~{t) < p~(s) < 1 mO 
for every t > s > . 

We next derive some compactness results for the sequences {w/c}fceN and {^fclfceN (defined by (|4.2| ). 

Lemma 4.11. Let {fc n }neN be the subsequence given by Lemma \4-10\ There exist a further subsequence (not 
relabeled) and a strongly measurable map u £ : [0, +oo) — > H 1 ^) such that Uk n (t) — u s (t) and Vk n {t) — ^ u £ (t) 
weakly in H 1 ^) for every t > 0. In addition, 

(i) u £ (0) = u a ; 

(ii) \\u £ {t)\\ L co(a) < *nax{\\uo\\L°°(0)> Il5l|z,°°(n)} for every t>0; 
(in) u £ g L°°(0, +oo; H\n)) H Lj^flO, +c^); H 2 (Q)) ; 

(iv) ^ = ml^O.+oo;^^)); 
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(v) u £ G AC 2 ([0, +00) ;L 2 (n)) and 

K(*)lli»(n) dt < ||Vuo||i2 ( n;R^) + ^ll^o - g\\h {Q) ; 

(vi) v' k — 1 u' e weakly in L 2 (0, +00; L 2 (f2)). 

Proof. We start by establishing the compactness of the sequence {ufc}. First Lemma [5751 yields for every i > 0, 

\\vk n (t)\\ L ~>(n) = ||wfcn(*)IU«(n) <max{||w |U~(n) ) ||fl , IU'»(n)}- (4-16) 

Then, combining the energy (in)equalities f|4. 1 1[) for Scheme 1 or (|4.12[) for Scheme 2 together with (|4.2p . we 
infer that 

sup ||Vufc„(t)|| L 2 (n;K ]v) < sup \\Vu kn (t)\\ L 2^ n . K N) < C £ , (4.17) 
t>o t>o 

for some constant C e > independent of fc„. Consequently, for every T > the set [J n Vk n ([0,T]) is relatively 
compact in L 2 (£l). On the other hand, (|4.11j) and (|4.12|l yield 

H n (r)\\ 2 mn)dr<£e(uo,p £ )<e s (uoA) < ||Vuo||£ a(n;HW) + ||K - fflll^n) . (4.18) 
Since for any t > s > we have 

lkW-»ik.Wll£»{n)< / H n ( r )hnn)dr, (4.19) 

^ s 

we deduce from (|4. 18[) and Cauchy-Schwarz inequality that 

WvkAt) ~ Vk n ( s )\\h(n) <{t-s) J \\v' kn (r)\\ 2 L2{n) dr < (t-s) ^||Vu ||!, 3(n . R w ) + ^\\u - g\\ 2 L2(n) 

By the Arzela-Ascoli Theorem we can find a subsequence of {k n } (not relabeled) such that 

v kn -> w e in ^°([0,T];L 2 (fi)) for every T > 0, (4.21) 

for some u e G < rf°' 1 ^ 2 ([0, +oo); L 2 (f2)). In particular, Ufe„(i) — > u E (t) strongly in L 2 (fl) for every t > 0, which 
yields (%) since Vk n (0) = Uo- 

On the other hand, in view of estimates (|4.16[) and ()4.17[) . we obtain (ii) and the fact that ffc„(i) — 1 u E (t) 
weakly in iJ 1 (f2) for every i > 0. By lower semicontinuity we also deduce from (|4.17j) that 

sup||« 6 (t)|| B i ( n) <C e . (4.22) 
t>o 

We next show the compactness of the sequence {uk n }- Let us now consider an arbitrary t > 0. For each neN 
there is a unique i G N such that t G (it jifc ]• We then have Ufc n (t) = Uk n {f k ) = Ufc„(i| )■ Consequently, by 

ll"fc„(i)-^(i)IU=(n) = H„(4j-u e (i)|U2 ( n) 

< K„(4„) ~ «fe„(i)IU 2 (o) + \\v kn (t) - u e (t)\\ L 2 {n} 



(4.20) 



<CvT*fcJ + IK(*)-«8(*)||i»(n) — >0. 

n— >oo 

Hence Uk n (i) — » u e (i) strongly in L 2 (O) for every i > 0, and in view of (|4. 1T[) we infer that Uk n (i) — *■ u e (i) weakly 
in i7 1 (il) for every i > 0. The mappings t h-> Uk n {t) being (locally) simple and measurable, we conclude as in the 
proof of Lemma 14.91 that u e : [0, +oo) — > 7? 1 (H) is strongly measurable. Moreover u e G L°°(0, +oo; 7J 1 (f2)) by 
(14.221) . By the pointwise strong L 2 (f2)-convergence of Uk n towards u e and the dominated convergence theorem, 
we have u kn —> u e strongly in L 2 OC (0, +oo; L 2 {Vt)). On the other hand, {uk n } is bounded in L 2 (0, T; H 2 {tt)) 
for every T > by (|4. 14[) . Hence itfc n — u £ weakly in L 2 OC (0, +oo; H 2 (£l)) which shows in particular that 
u e G L 2 OC (0, +oo; H 2 (£l)). Item (in) is thus proved. 

We now prove that u £ satisfies the Neumann boundary condition (iv). To this purpose let us fix T > and 
an arbitrary (p G L 2 (0, T; iJ x (0)). By Proposition we have = in L 2 (0, T; H 1 ^ 2 (dfl)), so that 

/ / (— Au E )ip dx dt = lim / / (— Auk )<pdxdt = lim / / Vit^ -Vipdxdt= / / Vu £ ■ V<p dx dt . 
Jo Jn »^°°Jo Jn n ^°°J Jn Jo Jn 
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From the arbitrariness of (p and T, we conclude that = in L 2 (0, T; H^ 2 (dfl)) for every T > 0. 

It remains to show the absolute continuity in time of u e . We note that (|4.18j) tells us that the functions 
Af~ n : t G (0, +oo ) i — ^ \\v' k (£)||i,2(q) are bounded in L 2 (0, +oo). Hence we can find a further subsequence (not 
relabeled) such that Ak n — 1 A weakly in L 2 (0, +oo), for a non-negative function A G L 2 (0, +oo) satisfying 

A 2 {t) dt < ||Vu || L 2 ( q. ]r n ) + -|K - g\\ 2 L 2 {n) . 

Letting n — > oo in (|4.19j) . we conclude that for every t > s > 0, 

\\u e (t) -u s (s)\\ L 2 {n) < J A{r)dr, 

which shows that u £ G ylC 2 ([0, +oo); L 2 (f2)), whence (v). 

Now, since is bounded in L 2 (0, +00; L 2 (f2)), up to a subsequence, } converges weakly in 

L 2 (0, +oo; L 2 (fl)) to some element in L 2 (Q, +oo; L 2 (f2)) which has to agree with it^. by (|4.21[) . This implies 
that (vi) holds. □ 

Remark 4.12. As a consequence of (Hi) and (v) in the previous lemma, u e : [0, +oo) — > is weakly 

continuous. 

We finally state a compactness result for the De Giorgi interpolants {u^j^gN (see Definition 14. lj) in the case 
of Scheme 1 . 

Lemma 4.13. If the discrete trajectories {(uk, Pk)}keN are obtained from Scheme 1, and {fc rl }„gN is the sub- 
sequence obtained in Lemma \4-ll\ then Uk n {t) u s (t) weakly in H 1 ^) for every t > 0. 

Proof. Let t > with t G ] • In view 01 63- (|47TT|) . and (j4~20]) . we have 



ll«ik„(<)-«fc„(*)IU»(n) < HfifcnW-Ufc^lU'cn) + IK-fc„W -^(C 1 )!^ 2 ^) ^ c v4^J, 

and thus Uk n (t) — > u e (t) strongly in L 2 (f2) by Lemma [4.111 On the other hand, (|4.7p and f|4. 1 1[) also show that 

II (i)||z,2(fJ;RJV ) < C £ , 

and thus Qfc„(t) — w e (t) weakly in i/ 1 (f2). □ 

As an immediate consequence of Lemmas 14.91 fc fl. Ill we obtain that unilateral (alternate or not) minimizing 
movements starting from (uq,Pq) do exist. 

Corollary 4.14. The collections GUMM(uq, p^) and GUAMM(uq,Pq) are not empty. 
5. Convergence of discrete trajectories 

The object of this section is to provide more accurate information on generalized unilateral (alternate) minimiz- 
ing movements, and to prove that they are solutions of the unilateral gradient flow of the Ambrosio-Tortorelli 
functional in sense of Definition 13.71 The main results of this section can be summarized in the following 
theorems. 



Theorem 5.1. Assume that f2 has a ff 1 ' 1 -boundary. Let uq G H 1 ^) n L°°(fi) and let p\ be given by (|3.1 
Any (u e , p £ ) G GUMM(uo, p%) is a curve of maximal unilateral slope for £ £ . More precisely, 

u e G AC 2 ([0, +oo); L 2 (0)) n L°°(0, +oo; H ^fi)) n if oc ([0, +oo); H 2 (il)) , 
p £ G L°°(0, +oo; W l ' p {il)) , < p E (t) < p E (s) < 1 for every t>s>0, 

and 

' u' e = Aiv({r, £ + p 2 £ )Vu £ ) -P{u e -g) in L 2 (0, +oo; L 2 (0)) , 

^ = in L 2 {0, +oo; H^dQ)) , 

u e (Q) = u , 
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with 

!£ £ (u £ (t), p e (t)) < £ £ (u £ (t), p) for every t > and p G W 1,v (Sl) such that p < p e (t) in ft , 
Pe(0)=p 6 O- 

Moreover, 1 1- > £ £ (u £ (t), p £ (t)) has finite pointwise variation in [0, +oo), and there exists an (at most) countable 
set jV e C (0, +oo) such that 

(i) (u E ,p e ) : [0, +oo) \7V £ — > x W 1,p (Vl) is strongly continuous; 

(ii) for every s G [0, +oo) \Af £ > a- n d every t > s, 

rt 

"" i £( r )lli2(n) 



£ £ (u £ (t),p £ (t)) + / \\u' e (r)\\ 2 mU )dr<£ e ( Ue (s),p e (s)). 



Theorem 5.2. Assume that has a < t? 1 ' 1 -boundary. Let uq G 1 (SI) nL°°(f2) and let p 1 ^ be given by (|3.1[) . Let 
(u £ ,p £ ) G GUAMM(uq,Pq), and assume that for each t > 0, (u £ (t), p £ (t)) is a strong L 2 (Q) x L p (Q)-limit of 
some discrete trajectories {(«&(£), Pk{t))}keN obtained from Scheme 2, and some sequence of partitions {^fc}/cgN 
of [0, +oo) satisfying \8k\ — > and 

( S ,+1 \ 

sup sup -|— < oo . (5.1) 

fceN \ i>i o k J 

Then all the conclusions of Theorem \5.1\ hold for (u e ,p e ). 



The entire section is devoted to the proof of those theorems. To this purpose, we consider for the rest of this 
section an open set f2 with 'rf 1,1 -boundary, an arbitrary element (u £ , p e ) in GUMM(uo, pff) or GU AMM(uq, pg), 
a sequence {<5fc}fc £ N of partitions of [0, +oo) satisfying \8k\ — > 0, and discrete trajectories {(«&, Pfe)}feeN associated 
to {<5fc}fc e N such that 

(u k {t),p k (t)) — > (u £ {t),p £ {t)) strongly in L 2 {tt) x D>(Q) for every t > 0. 

fc— >oo 

If (u e ,p £ ) G GUMM(uq,Pq) the discrete trajectories {(«&, Pfc)}feeN arc obtained from Scheme 1, while 
{(uk, Pk)}k£N are obtained from Scheme 2 if (u £ ,p £ ) G GUAMM(uq, pff). In addition, we require that (|5.ip 
holds if (u £l p £ ) G GUAMM(uq, pfj). Finally, extracting a subsequence if necessary, we may assume that all the 
results of Section |4] hold. For consistency we shall keep the notation of Section [U 

The plan of the proof is as follows. Let us first define the (diffuse) surface energy at a time t > by 

& s (t):= J(^\Vp s (t)f + ^(l-p e (t))A dx, (5.2) 

and the bulk energy 

*e{t):=\ f (Tfe + p 2 £ (t))Nu £ (t)\ 2 dx + | f (u £ (t) - gf dx . (5.3) 

In subsection 15. 1[ we establish a first minimality property for the phase field variable p £ leading to monotonicity 
and continuity properties of the surface energy & £ ( Proposition l5 . 3j ) . and then to the strong M^ 1,p (r2)-continuity 
of t H- p £ {t) outside a countable subset of (0,+oo). In subsection 15.21 we start by proving that u £ satisfies 
the inhomogeneous heat equation. Exploiting a semi-group property for this equation, we show that the bulk 
energy *8 e has also monotonicity and continuity properties (Proposition 15 -8[) from which the strong 7J 1 (i7)- 
continuity of t H> u £ (t) outside a countable subset of (0,+oo) follows (Corollary 15. 1Q[) . Then Subsection 15.31 is 
devoted to strong convergence results for the sequence {(uk(t), Pk(t))}k£^ in _ff 1 (i7) x M /1,P (J7) (Lemma 15.121 
& Proposition I5.13|) . As a consequence of these strong convergences, a stronger minimality property at every 
time for p £ is obtained. Finally, we show in Subsection 15.41 that (u £l p £ ) is a curve of maximal unilateral slope 
( Proposition 15 . 14)) and, as a byproduct, that an energy inequality is valid between (almost every) two arbitrary 
times (Corollarv l5.15j) . 
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5.1. Time continuity for the phase field variable and the surface energy 

We first establish several properties of the limiting phase field p e , starting from a (weak) minimality principle 
as stated in the following proposition. 

Proposition 5.3. For every t > 0, 

/ ( £ —\v Ps (t)\ p + ^-(i~Mt)) p ) dx<[ ffZ^\vpf + ^-(i- P y)dx 

Jn VP V £ / Jn V P P £ J 

for all p € W ' p (£l) such that p < p £ {t) in il. In particular, the surface energy & £ defined in (|5.2[) is non- 
decreasing on [0,+oo) 7 and thus continuous outside an (at most) countable set S £ C [0, +oo). 

Proof. Fix t > and let i E N be such that t £ (ijp 1 , 41- Consider a function p £ W 1,p (fi) such that p < p e (t) 
in f2, and define p k '■= p A Pk{t). Then pk £ W X ' P {Q) and pk < Pk(t) < p\ ■ By the minimality properties of 
the pair (v,k(t), pk(t)) obtained by either Scheme 1 or 2, 

£ e (u k (t),p k (t)) < £ £ (u k (t), pk) , 

and since p k < p k (t), 

I ( — \Vpk(t)\ p + ^(1 - p k (t))A dx< f ( — \Vpk\ p + 4(1 - Pk) p ) dx . (5.4) 
Jn \ P P £ / Jn \ P P £ / 

Let us now define the measurable sets Ak := {p < p k (t)}. By definition of p k , we have 

[\Vp k \ p dx=f \\7p\ p dx+ f \\7p k (t)\ p dx, 
Jn JA k Jn\A k 

and thanks to (|5.4[) . we infer that 

— / \Vp k (t)\ p dx+^- [ (1 - Pk (t)) p dx < — [ \Vp\ p dx+^- f(l- Pk ) p dx. (5.5) 
P JA k P £ Jn P Ja, P £ Jn 

Since p k (t) — > p £ (t) strongly in L p (fl) and p < p £ (t) in f2, we deduce that J£ N (Q. \ A k ) — > 0. As a consequence, 

f \Vp\ p dx^ [ \S7p\ p dx, 
J A k Jn 

and XA fc V ' p k {t) Vp e (t) weakly in L p (£l;M. N ) which in turn leads to 

liminf / \Vp k (t)\ p dx> [ \Vp £ {t)\ p dx. 
k ^°° JA k Jn 

Passing to the limit in (|5.5j) as k — > oo yields 

/ (—\Vp e (t)\ P + ^-(l-P S (t)) P ) dx< ( (— |Vp|* +^-{1-p)p) dx. 
Jn V P P e / Jn\ P P £ J 

In particular, taking p = p £ (s) with s > t leads to the desired monotonicity property for the function & £ . □ 

At this stage we do not have any a priori time-regularity for 1 1-> /? e (i) except that it is non-increasing, and 
thus it has finite pointwise variation (with values in L 1 (f2)). In the following result we show that this mapping 
is actually strongly continuous in M /1 ' p (0) outside a countable subset of (0, +oo) containing the discontinuity 
points of the surface energy & £ . 

Lemma 5.4. There exists an (at most) countable set TZ £ C (0,+oo) containing S £ such that the mapping 
t i — ^ p e (t) is strongly continuous in W /1,p (fi) on [0,+oo) \1Z £ . In particular, p £ is strongly continuous att = 0. 

Proof. Let TZ £ be the union of the set S £ given by Proposition 15.31 and the set of all discontinuity points of 

^ / p e (t)dx. (5.6) 
Jn 

Note that 1Z £ is at most countable by the decreasing property of the latter function. Let t £ [0, +00) \ TZ e , we 
claim that p £ is strongly continuous in W 1,p (il) at t. Consider a sequence t n —¥ t and extract a subsequence 
{t nj } C {t n } such that p £ {t n -) — 1 p± weakly in W 1,p (f2) for some p± £ W ' p (Cl). Upon extracting a further 
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subsequence, we may assume without loss of generality that t nj > t for each j € N (the other case t nj < t can 
be treated in a similar way). Then p e (t nj ) < p e (t) in f2, and passing to the limit yields p+ < p £ {t) in fl. On the 
other hand, by our choice of t as a continuity point of the mapping (|5.6[) . we have 

/ p e (t)dx= lim / p £ (t nj )dx = / p+ dx , 
Jn J^ifi Jn 

and thus = p e (t). As a consequence, the limit is independent of the choice of the subsequence, and the full 
sequence {p £ (t n )} weakly converges to p e (t) in W 1,p (n). Finally, using the fact that Ms a continuity point of & £ , 
we get that & £ (t n ) 6 £ (t), and thus \\p 6 (t„)\\w^p(n) -> llftC*)!!^.** (n)- We then deduce that p c (i„) p e (i) 
strongly in W 1>p (n). 

It now remains to show that p 6 is continuous at t = 0. Let t,- i be an arbitrary sequence. By Remark 14. 121 
we have u £ (tj) — 1 uo weakly in By Lemma WM p e G L°°(0, +00; M /1,p (0)), and we can extract a (not 

relabeled) subsequence such that p e (tj) — x p* weakly in W 1,p {Vt) for some p* G W 1,p (fi). According to the energy 
inequality proved in Lemma [4761 fc [4771 we have £ £ (uk{tj), Pkitj)) < £e(mo>Po) f° r au J G N and all k G N. Now 
we apply Lemma l3.1l to pass to the limit first as k ^ 00 and then as j — >• 00, which yields £ £ (uq, p*) < £ b (uq, Pq). 
From the minimality property (|3.1|) satisfied by , we deduce that £ £ (uq,P*) = £ £ (uo,Po)- By uniqueness of 
the solution of the minimization problem (|3.1[) . we have p* = Pq- Moreover, we infer from the discussion above 
that limj £ e (u e (tj), p e (tj)) = £ £ (uo,Pq), which implies that p £ (tj) —> P% strongly in W 1,p {Vi) by Lemma [3.11 
This convergence holds for the full sequence {tj} by uniqueness of the limit. □ 

Thanks to the just established continuity of t H- p s (i), we deduce that p~ and p e actually coincide almost 
everywhere in time whenever (|5.1|) holds. 



Corollary 5.5. Assume that (|5.1[) holds. Then there exists an J£ l -negligible set A4 £ C [0, +00) such that 
p £ {t) = p s (t) for every t G [0, +00) \ M £ . 

Proof. Step 1. Let us consider the function 4 : [0, +00) — > [0, +00) defined by 



'0 ifte[0,fj 



fej ! 



4(t) := < $i 

f * _1 + rf|r(*-4) ifte^t^withi^i 



Notice that 



Setting 



sup|4(*)-t| <3|* fc | — ► 0. 

t>0 fc^oo 



p*(t) := p £ (h(t)) , 

we infer from Lemma [5741 that — ► p £ (t) strongly in i 1 (f2) for every t € [0, +00) \ 7£ e . Since < p e < 1, by 
dominated convergence we have p\ — > p £ strongly in L 1 (0, T; L 1 (f2)) for every T > 0. Similarly, by p.5[) we have 
that pfc — > p £ and p,r — > p~ strongly in i 1 (0, T; i x (57)) for every T > 0. Given T > arbitrary, we estimate 

HftT - Pe||Li(0.T;Li(O)) < ||Ps - Pk II L 1 (0,T;£i (O)) + lift - (0,T;Z,i(n)) + llPe _ Af IUi(0,T;Li(n)) , > 0. 



Indeed, observing that p h (t) — pk(£k(t)) and £k{t) < £, we deduce from (|5.1D that 

lip* (<) - Pe(t)\\mu) dt= [ ||p fc (4(t)) - p e (4(*))IUi ( fi) A 

J51 




T 

\\pk(t)-p e (t)\\ LHa) dt—> 0. 



Hence |jp £ — p e ]|z,i(o, T ; z,i(m) = for every T > 0, whence p € (t) = p £ (t) for all t G [0, +oo) \ M E for some 
.^-negligible set A4 e C [0, +oo). □ 
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5.2. Time continuity for u £ and the bulk energy 

We start proving the following convergence result which is the main point to show that u £ solves the inhomo- 
geneous heat equation. 

Proposition 5.6. If(u e ,p e ) £ GUMM(u , pg) then 

div({r] e + pt)Vu k ) div((77 £ + p 2 £ )Vu £ ) weakly in £? oc ([0, +oo); L 2 (ft)) , 

while 

div({i h + (p-) 2 )\7u k ) div((?7 e + pf)Vu e ) wea% m L^ oc ([0,+oo);L 2 (fi)) , 

A;— >-oo 

i/( U£l p £ )6GMMM(ti 0) p§). 

Proof. 5*iep J. We first consider the case where (u £ ,p £ ) £ GUMM(uo, Po)- For t > such that t € (t^ -1 , t^] for 
some integer i > 1, the function Ufc(i) solves (|4. 15|) . Hence div((?7 e + p k (t))'Vu k (t)) £ L 2 (fi) with the estimate 

||div((% + pl(t))Vu k (t)) \\ L2{n} < \\v' k (t)\\ L2{n) + p\\u k {t) - g\\ L 2 m . 

In view of (|3.5|) and (|4.18|) . we deduce that for every T > 0, 

J \\dW((r h +pUt))Vu k (t))\\ 2 L2{n} dt<C J \W k {t)\\ 2 L2{n) dt + CT <C{T+l), 

for some constant C > independent of k and T. Hence we can find a subsequence (not relabeled) such that 
div((r? £ +p 2 k )\7u k ) ->> 6 weakly in L? oc ([0, +oo); L 2 (0)). Let ^ G ^ C °°([0,T] x fi). We have 

/ / div((r/ e + p 2 k )\Ju k )(pdxdt = - [ [ (n £ + p 2 k )Vu k -Vipdxdt . (5.7) 
Jo Jn Jo Jn 

From the convergences established in Lemmas 14.91 and 14.111 we have 

(»? e + p 2 k {t))Vu k {t) ■ V(p(t) dx — > (j} E + pl(t))Vu e {t) ■ V(p(t) dx for every t > 0. 



A;— f oo 

Then, from (|3.5j) and (|4.17|l we deduce that 



f [ (rj £ + p 2 k )Vu k -Vipdxdt — > [ [ + p 2 )\7u £ ■ \7tp dx dt , 
Jo Jn k^ccJ Q J n 



by dominated convergence. Therefore, letting k — > oo in (|5.7[) yields 

r T r i-T 



f I Qtpdxdt= - ( ( (r) e + p 2 £ ) Vu £ ■ Vip dx dt 
Jo Jn Jo Jn 



from which we deduce that O = div((?7 £ + p 2 )VM e ), and this first step is complete. 

Step 2. The case where (u £ ,p £ ) € GUAMM(uo,p^) is essentially the same, replacing p k by the translated p k , 
and using the property established in Corollary 15.51 We do not reproduce the proof for this case. □ 

Corollary 5.7. The function u £ G AC 2 ([0, +oo); L 2 (il)) solves 

' u' £ = div((r/ e + p 2 £ )Vu s ) - f3(u £ - g) in L 2 (0, +oo; L 2 (n)), 

5^ = mL^O.+oo;^^)), (5.8) 

ov 

u £ (0) = u Q . 

Proof. It suffices to combine (|4.15[) (replacing p k (t) by p k {t) for the case (u £ ,p £ ) £ GU AMM(uq, Pq)) with 
the convergence results of Lemma [4.111 and Proposition l5.6l □ 
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We are now in position to prove the decrease nionotonicity of the bulk energy *B e . 



Proposition 5.8. Let to > and set p £ °{t) := p E (t + to). For any wq G fl L°°(tt) there exists a unique 

solution w £ G AC 2 {[0, +oo); L 2 (0)) D L°°(0, +oo; fl^fi)) of 



'w' £ = div((r ?e + (p*°) 2 )Vu; £ ) - /3(w e -5) L 2 OC ([0, +00); H' 1 ^)) , 
< ( Ve + (p £ °) 2 )Vw £ ■ 1/ = in Lf oc ([0, +oo); ff-Va^fi)) , 

and iu e satisfies the following energy inequality for every t > 0, 



(5.9) 



(r l£ + (p t £ °(t)) 2 )\Vw £ (t)\ 2 dx+t- / {w E {t)-gfdx< - / (7 ?£ +p 2 (io))|V™ | 2 / (u; -ff) 2 ^. (5.10) 



/3 



In particular, for any t > i/ie function u E { ■ +to) is the unique solution of (|5.9[) w&h initial datum wo '■= u £ (to). 
As a consequence, the bulk energy Q5 e defined in (|5.3p zs non increasing on [0, +00), and £/iits continuous outside 
an (at most) countable subset B e of [0, +00). 

Proof. Step 1, Uniqueness. Let w £t i and w £; 2 be two solutions of (|5.9p . and set z e := w £ ^\—w £> 2- Then z £ (0) = 0. 
The variational formulation of (|5.9p implies that for any T > and any test function <j> G L 2 (0, T; iJ 1 (SI)), 

/ / (z' e 4>+ (rj E + (pl° ) 2 )Vz £ • V0 + I3z £ <t>) dxdt = Q. 
Jo Jn 

Choosing <j>(t) := z £ (t)x\o,T](t) as test function above yields 

/ [ z' e z £ dxdt < for every T > . 
Jo JO 



On the other hand, since z £ <E AC 2 ([0, +00); L 2 (Q)), we have ||^ e (-)l| 2 im) 6 AC ([0, +00)) and 

d 



Therefore, 



> 



||ze(<)||£2(n) = 2 / z' £ (t)z £ (t)dx for a.c. t £ (0, +00) 



z^z £ didt = — ||z E (T)||jr2m') for every T > , 



'0 Jo 

which shows that z £ = 0, i.e., w Et i = w £ _2- 

Step 2, Existence. For what concerns existence, we reproduce a minimizing movement scheme as before. More 
precisely, given a sequence ! 0, we set := irfe for i € N. Taking w® := tooj we define recursively for all 
integer i > 1, w k £ H (Q) as the unique solution of the minimization problem 

n e + (p^(tr 1 )) 2 )|Vz;| 2 ^ + ^ / + J- / {v-w^fdx 



mm 

•jGH 1 (n) [2 



k Jn 



Using the minimality of w k at each step and the fact that < /Q*°(r^) < p £ °{r k 1 ), we obtain that for every 
integer i > 1, 

J / ( ?? e + ^°(rr 1 )) 2 )|V^| 2 d a; + ^ / ( w * - 5 ) 2 d.T + ]T-i- /" (njl-^fdx 



1 



<7^ / (V, + p;(to))\Vw \' dx + ^ / (t«o-flrdaj. (5.11) 

Let us now define the following piecewise constant and piecewise affine interpolations. Set Wfc(0) = Wk(0) = Wq, 
and for t £ (t£~\t|], 



w k (t) 

Qi°(t) 

u) fc (t) 



= P^(^- 1 ), 



< 1 +r i r 1 (t-rj- 1 )K-< 1 ) 
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By Lemma I5T41 we have Q k °(t) —¥ p £ °(t) strongly in W 1,p (fl) for all t G [0, +oo) \ (— 1 + 1Z £ ). Arguing exactly as 
in the proof of Lemma 14.111 we prove that (for a suitable subsequence) Wk(t) w £ (t) weakly in H 1 ^) for every 
t > and w' k w' E weakly in L 2 (0, +oo; L 2 {Q)), for some w £ G AC* 2 ([0, +oo); L 2 {Q)) n L°°(0, +oo; ff 1 ^)) n 
^1oc([0j +°°); -ff 2 (^))- Then we can reproduce with minor modiheations the proof of Proposition HT61 and Corol- 
lary 15.71 to show that w £ is a solution of (j5.9|) . 

Since < p*°(t) < Qu{t) and Wfc(t) — > w s (t) strongly in L 2 (Q) for every t > 0, we infer from (|5.11l) that for 
every t > 0, 

\ [ (V* + p 2 e(to))\Vw Q \ 2 dx + | f {w -g) 2 dx 



> lim 



iminf Q ^ + (fi k ° (t)) 2 ) \Vw k (t)\ 2 dx + | jf M*) - ff) 2 

>hminfi / ( ??£ + (p*"(i)) 2 )|V^W| 2 ^ + ^ / (« £ (t)-5) 2 ^ 
fc->oo 2 J n L J n 

>\f {vs + (p £ a (t)) 2 )\Vw £ (t)\ 2 dx+^ f (w e (t) - g) 2 dx , 

and (|5.10j) is proved. □ 

Remark 5.9. We notice that the proof of Lemma [5^41 together with Remark 14.121 show that the function ( B £ is 
actually continuous at time t = 0, i.e., g" B £ . 

As a consequence of Lcmma l5.4l and Proposition l5.8[ we obtain the strong continuity in iJ 1 (f^) of the mapping 
t t— > u £ (t) outside a countable subset of (0, +00) containing the discontinuity points of & £ and ( 3 £ . 

Corollary 5.10. The mapping u £ : [0, +00) — > is strongly continuous on [0, +00) \ (1Z £ U B £ ). 

Proof. Let us consider to G [0, +00) \ (1Z £ U B £ ) and {t„} C [0, +00) an arbitrary sequence such that 
t n to- Since t a ^ U £ U B £ we have < B £ {t n ) -> <8 e (io) and p £ (t„) -> /0 e (*o) strongly in W 1,p (fi). Therefore 
£ £ (u £ (t n ), p £ (t n )) — > £ E (u e (to), p e (to))- On the other hand u £ (t n ) — u e (io) weakly in by Remark [4. 12[ 

and the conclusion follows from Lemma T3. II □ 

5.3. Strong convergences and limiting minimality 

Thanks to the equation solved by u e , we are now able to improve the weak i? 1 (r2)-convergence of the sequence 
{wfcj-fcgN hito a strong convergence. We start by proving that the bulk energy converges in time average. 

Lemma 5.11. If (u £ ,p £ ) G GUMM(uo, Po) we have 

lim / / (r) £ + p 2 k (r))\S7u k {r)\ 2 dx dr = [ f (r) £ + p 2 £ (r))\Vu £ (r)\ 2 dx dr for every t > s > , (5.12) 



while 

lim / / (r/ £ + (p k (r)) 2 )\\7u k (r)\ 2 dx dr = / / {r\ £ + p 2 {r))\Vu £ {r)\ 2 dx dr for every t > s > , 
k^^Js Jn Js Jn 

tf(u £ ,p £ )eGUAMM(u ,p £ o). 

Proof. We only consider the case (u £ , p £ ) £ GUMM(u , Po)- In the case (u £ , p £ ) G GU AMM(uo, pf,), it suffices 
to reproduce the argument below with p k instead of p k . 

Taking u k (r) as test function in the variational formulation (|4.15|) and integrating in time between s and t 
leads to 

(r] £ + p 2 k (r))\Vu k (r)\ 2 dxdr = - / / v' k (r)u k (r) dx dr - /3 / / (u k (r) - g)u k (r) dx dr . 
n Js Jn Js Jn 

From Lemma 14.111 we have u k — > u £ strongly in L 2 OC ([0, +00); L 2 (il)) and v' k — ^ u' £ weakly in L 2 (0, +00; L 2 (Sl)). 
Therefore, 

lim / / (r) e + p 2 k (r))\Vu k (r)\ 2 dx dt = - / / u' £ (r)u £ (r) dx dr ~ (3 / / (u £ (r) - g)u £ {r) dx dr . 

Js Jn Js Jn 
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On the other hand, according to equation (|5.8[) solved by u e , we have 

u' e (r)u e (r) dx dr — fi I I (u e (r) — g)u e (r) dx dr = / / (r] e + p 2 (r))\Vu E (r)\ 2 dx dr , 
s Jn Js Jn Js Jn 

which leads to (|5.12[) . □ 

From Lemma l5.11[ we deduce the announced strong convergence of the sequence {uk}k£N- 
Lemma 5.12. For every t € [0, +oo) \ B e , u k (t) — > u £ (t) strongly in if 1 (SI) . 

Proof. Step 1. We first assume that (u e ,p e ) G GUMM(uq, py). Let to S [0, +oo) \ B e . Since Q3 6 is continuous 
at t , for every a > there exists S a > such that Q5 e (i) < 25 £ (io) + a for all t G [to - S a ,to\- 

Let us fix a > arbitrary. Since £ £ (u k ,p k ) < E^u 1 ^ 1 , p k ) and p l k < p\~ x in fl for each integers /c and i > 1, 
we infer that the function 

f ^\l (Ve+pl(t))\Vu k (t)\ 2 dx+^ f (u k (t)-g) 2 dx 
is non increasing on [0, +00), and thus 

Sa Q jjjle + p 2 k {t ))\Vu k {t )\ 2 dx + ^ jf (ukfo) - 3) 2 da: 

< f J / (% + p 2 k (t))\Vu k (t)\ 2 cfe + § / (u k (t) - 

Jto-s a \ z Jn L Jn 

By Lemma \5. Ill and the strong convergence of u k to u £ in Li 2 oc ([0, +00); L 2 (il)), we infer that 

5 Q limsup (i I + p 2 k (t ))\\7u k (t )\ 2 dx + ^- [ (u k (t a ) - g) 2 dx 
fc->oo \^ Jn ^ Jn 

< lim f (\( + ^(t))| V Wfe (t)| 2 dx+^- f (u k (t) - g) 2 dx) dt 
k ^°° Jt -s a \ z Jn 1 Jn J 

to 

"& e {t)dt < (m E (t )+a)S a . 

to-Sa, 

Dividing the previous inequality by S a and using the strong convergence of u k (to) in L 2 (fl), wc derive in view 
of the arbitrariness of a that 

limsup / (r le + p 2 k {t Q ))\\7u k (t )\ 2 dx < / ( Ve + p 2 {t a ))\\7u £ {t )\ 2 dx . 
fe->oo Jn Jn 

As in the proof of Lemma 13.11 we obtain 



(r le + pi(t ))\Vu £ {t )\ z dx <limM / (r? £ + pl(t ))\Vu k {t )\ 2 dx . 

fe ^°° Jn 

Combining the last two inequalities we conclude 

lim / (r/ £ + p 2 k (tn))\\7u k {t )\ 2 dx = [ {r h + p 2 £ (t ))\\7u e {t )\ 2 dx . (5.13) 
fc ^°° Jn Jn 

Finally, using (|5.13[) , the weak convergence of p k {to) to p 6 (t ) in M /1 ' P (J7), and the weak convergence of u k (to) 
to u e (to) in H 1 ^!), we can argue as in the proof of Lemma [37X1 to show that u k (tn) —> u £ (tn) strongly in iJ 1 (f2). 

Step 2. If (u E ,p e ) e GUAMM(uq,Po), we essentially proceed as in Step 1 up to the modification indicated 
below. First notice that inequality (|4.13|) shows that for each integer k the function 

\ I (r ] e + (p k (t)) 2 )\Vu k (t)\ 2 dx + ^ I (u k (t)-g) 2 dx 

is non increasing on [0, +00). Then we conclude as in Step 1 that 

limsup / (r? £ + p 2 k {t ))\Vu k {t Q )\ 2 dx < limsup / (r, s + (p k (t )) 2 )\Vu k (t )\ 2 dx 

fc->oo Jn fc->oo 



< / (n £ + p 2 e (t ))\Vu £ (t )\ 2 dx . 



30 J.-F. Babadjian & V. Millot 



At this stage, it suffices to continue the argument of Step 1 to reach the conclusion. □ 

We now derive the (strong) minimality property for p e (t) at all times, as well as the strong W 1,P {£1)- 
convergence of {p k {t)}k&i at all continuity points of the bulk energy. 

Proposition 5.13. For every t > the function p e (t) satisfies 

£ e (u £ {t), p e (t)) < S s (u s (t), p) for all p G W 1 ' P (Q) such that p < p e (t) in Ct . (5.14) 

In addition, if t € [0, +oo) \ B e then 

p e (t) = argmin {£ e (u e (i),p) : p G W 1 ^^) such that p < p~(t) in tt} , (5.15) 

and Pk(t) — > p e {t) strongly in W 1,P (Q). 

Proof. Let us fix an arbitrary t > 0. Since u k {t) — ^ u e (t) weakly in iJ 1 (f2), we can find a (not relabeled) 
subsequence and a nonnegative Radon measure p G ^#(M W ) supported in Vl such that 

\Vu k (t)\ 2 ^ N L ft -» \\7u s {t)\ 2 ^ N \_n + p 

weakly* in jtf(R N ). Then we consider the functionals Tk and T defined on W 1,p (£i) by 



, = :£ £ (u k (t),p) if p < Pk -(t) , = I £e(Mt),P) + \ IJVs + (?) dp if p < Pe ( 



otherwise , 

I +co otherwise . 

Since pl(t) p £ (t) weakly in W 1,p (Vl), we may argue as in the proofs of Proposition 13.91 (Step 1) and 
Proposition 13. 1 II to show that Tk T-converges to J- with respect to the sequential weak topology in W 1,p (il). 
Since 

p k (t) = argmin F k (p) , 

and pk(t) — 1 Pe(t) weakly in W 1,P (Q), we infer that 

p e (t) = argmin J-(p) . (5.16) 

Let us now fix an arbitrary p G W 1,p {Vt) such that p < p e (t) in f2, and set p + := p A 0. Then p + G W 1,p (f2), 
< p + < p E {t) in Q,, and £ e (ii £ (i), p+) < £ £ {u E {t),p). Since p+ < p £ (i) < p~(t) in ft, we have T{p E {t)) < T(p+) 
which leads to 

S s (u £ (t),p e (t)) < £ e (ue(t), Pe (t)) + \j {p 2 M - (p + ) 2 )dp < £ e (u e (t),p+) < £ e (u £ (t),p), 
and (|5.14[) is proved. 

Next we observe that if t G [0,+og) \ B e , then p = by Lemma [5.121 Hence J-(p) = £ £ {u e (t),p) for every 
p G W l ' p (Q) such that p < p7(t) in 0, and (|5.15[) is a consequence of (|5.16[) . From the F-convergence of F k to 
J- we also have min J- k —> min J 7 , and thus 

£ s {u k {t), Pk {t)) — > ^(« f (i), p e (t)) , 
and the strong convergence in VF 1,p (f2) of Pk(t) follows from Lemma 1331 □ 

5.4. Energy inequalities 

Proposition 5.14. The mapping (u £ ,p e ) : [0, +oo) — > L 2 (fl) x L P (Q) is a curve of maximal unilateral slope 
for £ e . 

Proof. Let us define for each k G N and t > 0, \ k (t) := £ £ (uk(t), pk(t)). By Lemmas 14.61 fc 14771 the function 
Afc : [0, +oo) — [0, +oo) is non-increasing and bounded uniformly with respect to k. By Helly's Theorem for 
monotone functions we can find a (not relabeled) subsequence such that 

Afc(i) — > \(t) for every t > , 
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for some non-increasing function A : [0, +00) — > [0, +00). Then we infer from Lemma 15.121 and Proposition 15 . 13l 
that 

X(t)=£ e (u e (t),p c (t)) for every t G [0, +00) \ B £ . (5.17) 

We shall now distinguish both cases (u £ ,p £ ) G GUMM(uq, p$) and (u e ,p e ) £ GUAMM(uo,Po). 

Case 1. We first assume that (u e ,p e ) G GUMM(uo, Pq). Let us fix t > s > 0, and write s £ {t\ s ~ ,tl s ] and 
i G (^fe _1 , ife ]■ As in the proof of Lemma B~B1 we deduce from f|4.8[) that 

Afc(*) + J f" H(r)\\% m dr+ \ f " \G k (r)\ 2 dr < \ k (s) . (5.18) 



Applying (vi) in Lemma [4.111 and Fatou's lemma, we let k — > 00 in (|5.18[) to obtain 

m + \f s IKWII' 3 («) dr + 5 / fe** ' Gfc Wl 2 dr ^ A ( fi ) ■ 

We now claim that there exists an Jz? 1 -negligible set C £ C (0, +00) such that 

liminf \G k {t)\ 2 > \d£ £ \ 2 (u £ (t),p £ (t)) for every t G (0,+oo)\£ e . (5.19) 

fe— >-CX3 

Before proving this claim we complete the argument. Thanks to (|5.19p we have 

A(i) + ^ IKMII^(n) dr + \ J \d£ £ \ 2 (u £ (r),p £ (r)) dr < X(s) . (5.20) 

Notice that the second integral term above is well defined by the explicit formula (|3.10p . Proposition [5J21 and 
items (Hi) & (iv) in Lemma T4. Ill Since s and £ are arbitrary, we deduce from (|5.20p that 

A'(f) < ~H(t)\\ 2 L 2 (Q) - i|a£ £ | 2 ( U£ (i),PeW) for a.c. t G (0,+oo) . 

In view of (|5.17|) . we have thus proved that (u £ ,p £ ) is a curve of maximal unilateral slope for £ £ . 
It now remains to prove (|5 . 19[) . First notice that Lemma |4"U1 implies through Fatou's lemma that 

r+oo 

/ liminf \G k {r)\ 2 dr < +00 . 
Jo fc->oo 

Hence we can find an JC -'-negligible set L £ C (0, +00) such that 

C(t) := liminf \G k {t)\ 2 < +00 for every t G (0,+oo)\£ e . 

We set C £ := C £ U B E , so that ^f 1 (£ e ) = 0. Let us now fix an arbitrary t G (0, +00) \ C £ . We then extract a 
subsequence {fc„} (depending on t) such that 

lim |G fc „(i)| 2 = C{t). 

n— yoo 

In view of estimate (|4.9p we deduce that 

sup|9£ e |(%„(i),p fc „(t)) < +00. (5.21) 

new 

On the other hand, we infer from (|4.7p and (|4.1ip that 

sup £ju kn (t),p kn (t)) <+oo. (5.22) 

rigN 

Then, (|5.2ip together with (|5.22p implies that the sequence {u kn (t)} is bounded in H 2 (il) by Proposition 13.91 
and thus U kn {t) —¥ u £ (t) strongly in i? x (0) by Lemma |3~T51 By ([5T2"2"j) . {p kn (t)} is bounded in and 
up to the extraction of a further subsequence we may assume that p kn {t) — v p* weakly in M /1,p (f2) for some 
p* £ W^ifl). 

Now we consider once again the functionals T k and J- defined on W 1,p (f2) by 

«^*w, Md ^ )! _( 4( "- W '" ,f ^"- (i) - 

I+ 00 otherwise, +00 otherwise. 
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Arguing as in the proof of Proposition 15.131 we obtain that J- k T-convcrges as k — > oo to J- with respect to the 
sequential weak topology of W 1,p (d). By the very definition of Pk n {t) we have 

pk n (*) = argmin T kn (p) , 

pew l .p(n) 

and since Pk n {t) — ^ P* weakly in W 1,p (f2), we deduce that 

p* = argmin F(p) . 

P ew 1 .p(Q.) 

Since t g" B £ we conclude from Proposition 15.131 that p* = p £ (t) by uniqueness of the minimizer (again due to 
the strict convexity of £ £ (u £ (t), •)). 

To summarize, we have established that Uk n {t) — > u £ (t) strongly in and that pfc n (t) — Pe(t) weakly 

in W 1,p (Vl). By (|5.21[) we can now apply Proposition 13.111 and use (|4.9[) to derive that 

\d£ £ \ 2 {u £ (t),p £ {t)) <liminf \d£ £ \ 2 {u kn {t),~Pk n (t)) < Km \G kn {t)\\ 

n— too n— too 

and (|5 . 19[) is proved. 

Case 2. We now assume that (u £ , p £ ) £ GUAMAI(u , pg). Let us fix again t > s > 0, and write s G (ii" -1 ,^] 
and t G (ifc 1 !^*]- From Lemma T4. 71 we deduce that 

A fe (t)+/ ll<M|| 2 L2( n)^< A fc (s). (5.23) 

Still applying in Lemma [4.111 and Fatou's lemma, we let k —> oo in (|5 . 23[) to obtain 

A(t)+ / ||<(r)||| 2(n) dr<A( S ). (5.24) 

On the other hand, in view of Lemma 14.111 (iii)-(iv), the minimality (|5 . 14[) satisfied by p e , and the equation 
(|5.8[) satisfied by u e , we infer from Proposition 13.91 that 

H(r)\\ L , m = \d£ £ \(u £ (r),p e (r)) for a.e. r e (0, +oo) . (5.25) 

Then (|5.24p yields (|5.20[) , and the conclusion follows as in Case 1. □ 

Corollary 5.15. For every s £ [0, +oo) \ B £ and every t > s, 

£ £ (u e (t), p £ (t)) + J K(r)||£ a(n) dr < £ £ (u e (s),p e (s)) . 

Proof. Step 1. We assume in this step that (u £ ,p £ ) £ GUMM(uq,Pq). Let us fix s e (0, +oo) \ B £ and t > s. 
Using the notation of the previous proof, we first notice that £ £ (u £ (t), p £ {t)) < X(t) for every t > by Lemma [3~T1 
Then, combining (|5~T7)) with (|5T20|) yields 

£ £ (u £ (t),p £ (t)) + X - J \K(r)\\ 2 L2{n} dr+\J^ \d£ £ \ 2 (u £ (r), p £ (r)) dr < £ £ (u £ (s), p £ (s)) . 

As in the previous proof, Lemma \4.11\ (iii)-(iv). ()5.14p . (|5.8|) . and Proposition 13.91 imply that (|5.25j) holds, and 
the result follows. 

Step 2. If (u e ,p e ) G GUAMM(u ,pl), we still have £ £ {u £ (t), p £ {t)) < X(t) for every t > by LcmmaO Then 
the result is a direct consequence of (|5.17[) and (|5.24j) . □ 



6. Asymptotic behavior of unilateral minimizing movements in the Mumford-Shah limit 

The main goal of this section is to analyse the behavior of a unilateral (alternate) minimizing movement as 
£ tends to zero. We prove that in the limit e — > 0, we recover a parabolic type evolution for the Mumford- 
Shah functional under the irreversible growth constraint on the crack set similar to [12] . The result rests on the 
approximation of the Mumford-Shah functional by the Ambrosio-Tortorclli functional by means of r-convergence 
proved in |5|6I28) . The main result of this section is the following theorem. 

Theorem 6.1. Assume that Q has a ff 1,1 -boundary. Let e n 1 be an arbitrary sequence, uq G H (CI) r\L°°(p,), 
and Pq™ determined by (|3.1[) . Let {(u £n , p £n )} n£ jq be a sequence in either GU MM (uq, pg") or GU AMM(uq, Pg"). 
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In the case {(u Sn , Pe„)}neN C GUAMM(uo,pQ n ), assume in addition that assumption (|5.1[) in Theorem \ 5.2\ 
holds. Then there exist a (not relabeled) subsequence and u* G ^4C 2 ([0, +00); L 2 (Cl)) such that 

Pen(0 — ^ 1 strongly in L p (Cl) for every t > , 

u e n {t) ~~ ^ u*{t) strongly in L 2 (Cl) for every t > , (6-1) 
u' £ v!^ weakly in L 2 (0, +00; L 2 (Cl)) . 

For every t > the function u*(t) belongs to SBV 2 (Cl) n L°°(Ct) with 

IK(*)IU°°(fi) < max {|l M ollL-(o), IblU^cn)} , (6-2) 
andVu* E L°°(0, +00; L 2 (£l; R^)). Moreover?/* sofoes 

{< = div(Vu*) - - g) in L 2 {0, +00; L 2 {Cl)) , 
Vit*-^ = mL 2 (0,+oo;i/'- 1 / 2 (an)), 
U»(0) = M , 

and there exists a family of countably J$? N ~ 1 -rectifiable subsets {r(i)} t >o of CI such that 

(i) r(s) C r(i) for every < s <t; 
(H) Ju,(t) C r(t) /or every t > 0; 
(Hi) for every t > 0, 

i / |V^(i)| 2 d 2; + ^ Ar - 1 (r(t)) + ^ f (u4t)~g) 2 dx+ [ H(s)\\ 2 L2{n) ds 

|Vn | 2 d;r + ^ / (u - 5 ) 2 da;. 

This section is now essentially devoted to the proof of this theorem. To this purpose, we consider for the rest 
of the section an open set CI with ^ 1,1 -boundary, a sequence e„ 1 0, and an arbitrary sequence {{u £n , /0 e „)}neN 
in GUMM(uo,Pq) or GU AMM{u$, p§) (and we assume that the assumptions of Theorem 15.21 hold in this later 
case). 



6.1. Compactness and the limiting heat equation 

We start by proving compactness properties for the sequence {(u £n , p £ll )} ng N- 

Proposition 6.2. There exist a (not relabeled) subsequence {u ere }„ 6 N and a function u* G AC 2 ([0, +00); L 2 {Cl)) 
such that dHH) holds. In addition, u*(t) G SBV 2 {Cl) n L°°(0) with $£2$ for every t > 0, and i/ie mapping 
t h-> Vu*(i) G Z/^fijR^) is strongly measurable with Vu* G L°°(0, +00; L 2 (S1; E^)). Moreover, /or every t > 
and any < <5i < $2 < 1, there exists s n = s n {t, 61,62) G (^1,^2) swc/i £/iai £?„ := {p £n (t) < s "} ^ as finite 
perimeter in Cl, u £n {t) := (1 - XE n )u Sn (t) G SBV 2 {Ct) n L°°(Cl), and 

Ue n (t) ~~ ^ w*(i) strongly in L 2 (Cl) , 
u £n (t) — weakly* in L°°(Cl) , 
Wu Sn (t) Vu*(i) wea£% ^(fijK^) . 

Finally, for any open subset A C Cl, 

< , . A (6.3) 

/ |V«»(t)| 2 dx < liminf / (n e + p 2 (i))|Vu e (£)| 2 da; . 
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Proof. Step 1. We first derive a priori estimates from the energy inequality obtained in Corollarv l5.15l Indeed 
according to that result together with the minimality property (|3.ip of pg™, we infer that for every t > 0, 

1 / (t))\Vu en (t)\ 2 dx+^ { \V Pen (t)\*dx+-2- [ {l- Pen {t)Ydx 



P Jn n P'£n Jn 

+ [ H n (s)f L2(Q) ds<£ en (u ,p e n )<S en (u ,l)< I \Vu \ 2 dx + ^ [ (u Q -g) 2 dx. (6.4) 
Jo Jn 1 Jn 

Then, applying Young's inequality and using (|2.2[) . we obtain 

^ / \Vp Sn {t)?dx + ^- { {l-p £n {t)fdx> P - f (l-p En (t)y-'\Vp En (t)\dx, (6.5) 
P Jn P £ n Jn 2 J n 

from which we deduce the following uniform bound 

IK lliwoo^n)) + f(l°~ + Pi W)l Vu *» Wl 2 dx + JO- - (0) p_1 |VPe„ (t)\dx < Co , (6.6) 
for some constant Co > independent of e n and t. 

Step 2. We now establish the weak convergence of {u £rl } and the bound (|6.2[) . Recalling that < p Sn < 1, the 
fact that 

p s {t) — > 1 strongly in L p (fl) for every t > , 

is a direct consequence of (|6.4[) . According to items (v) and (m,) in Lemma 14.111 the sequence {u £n } in uniformly 
equi-continuous in L 2 (17), and for each t € [0, +oo), the sequence {u £n (t)} is sequentially weakly relatively 
compact in L 2 (fl). Therefore, according to Ascoli-Arzela Theorem, we can find a (not relabeled) subsequence 
and S AC 2 ([0, +co); L 2 (fl)) such that u £n (t) — 1 u* (i) weakly in L 2 (f2) (and also weakly* in L°°(Q)) for every 
t > 0, and — ^ «'„ weakly in L 2 (0, +oo; L 2 (f2)). In particular, (|6.2p follows from item (zi,) in Lemma [4.111 

S^ep 5. We now examine more accurately the asymptotic behavior of the sequence {u £n } as in [516128] . and 
prove (|6.3[) . Let us fix t > 0, < Si < 82 < 1 and an arbitrary open subset A of f2. According to the £?V^-coarea 
formula (see [3j Theorem 3.40]), 

(1 - p £n (i)f- 1 \Vp £n (t)\ dx = [ (1 - s)P- 1 Jf N ~ 1 (d*{p £n (t) < S }f)A) ds 

Jo 

> f \l-s) p - 1 Jif N - 1 (d*{p en (t) <s}nA)ds. (6.7) 
J s t 

Consequently, by the mean value theorem there exists some s n = s n (t, 81, 82, A) £ {8\, 82) such that 

8 P - 8 P 

(1 - Pe n (t)) p ~ 1 \V p £n (t) \ dx > — — —Jif N ~ 1 (d*E n n A) , (6.8) 

P 

where E n := {p £n (t) < s n } fl A. Note that from (|6.4[) we have 

JZ N (E n )<- / (1 - p £n (t)T dx < " ^0asn->oo, (6.9) 

for some constant C > independent of n. 
Let us define the new sequence 

u £n {t) := {I - X E n )u £n {t) ■ (6.10) 

By (|6.9[) we have 

IKJO-u^lU^) < || Ue j L » (n) ^/.s?"(£ n ) ->o (6.11) 

as n — > 00, from which we deduce that u £n (t) — ti*(i) weakly in L 2 (A). On the other hand, according to [3J 
Theorem 3.84] we have u £n (t) G SW 2 (A) n L°°(A) with 

Jfi en (t) C C**-E n , 

Vfi en (t) = (l-XB»)Vu B „(t) ) 
Jwe„(*)IU«>(A) < max{||uo|| £«>(«), ||ff||i=o (n )} . 
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By the energy estimate (|6.6[) together with (|6.7p and (16.81) . 

llvu Sn (t)ni 2(A . KN) < i / ^(^ivu^^i 2 ^ < § 



and 



^- 1 {Ju en( t ) nA)<^- 1 (d*E n nA) 



< 



CoP 



We are now in position to apply Ambrosio's compactness Theorem in SBV (see Theorems 4.7 and 4.8 in [3]) 
to deduce that u*(t) g 5By 2 (r2) (by arbitrariness of A), and that 

u £n {t) — >• u*(i) strongly in L 2 (A) , 
u E „(i) weakly* in L°°(A) , 

Vu e „(t) Vu*(t) weakly L 2 (A; R w ) . 

In view of (|6. 1 1 1) we deduce that u £n (t) —> u*(t) strongly in L 2 (n) for each t > (again by arbitrariness of A). 
Next Proposition 12. II yields 

2M ,N - l {J u u ) nA)< \immiJ^ N - 1 (d*E n n A) . 

n— >oo 

Combining this inequality with (|6.8[) wc get that 

(^-^)jr JV - 1 (J„. (t) nA)<limiiif I / (l-p^W^IVp^Wldx, 

Tl->OC 2 J A 

and the first inequality of (|6.3[) follows by letting <5i — >■ and 62 — > 1. 
For what concerns the bulk energy, we have 

{i] £n + p 2 n (t))\\7u £n (t)\ 2 dx > s 2 n \Vu £n (t)\ 2 dx > 5\ I \Vu £n (t)\ 2 dx . 

J A\E n J A\E n 

Since u £n (t) = (1 — xb„)w £?i (<), we have Vu £n (t) = (1 — XE n )Vue n (t) Vu*(t) weakly in L 2 (A;M. N ), and thus 

liminf / (r) E + p 2 (i))|Vu e {t)\ 2 dx > 5\ liminf / \\7u £ {t)\ 2 dx>8l [ |Vu»(i)| 2 dx , 
and the second inequality of (|6.3p follows by letting <5i — >• 1. 



Step 4- It now remains to prove the strong measurability in L 2 (il;M. N ) of t >-> Vu»(t), and that Vu» € 
L°°(0, +00; L 2 (S7; R N )). Given t > and < <5i < 52 < 1 arbitrary, let us consider as in Step 3 the set 
E n and the function u Sn (t) <E SBV 2 (n) given by (jOO]) with A = Cl. Then, 

(%„ +P?„ (*))(! -XB„)Vu eB (t) = (lfe„ +pg n (t))Vfi Bn (t). 

Note that this last sequence is bounded in L 2 (f2;IR Ar ). Since p £n (t) — ► 1 strongly in L p (tt) with < p £n < 1, 
and Vu £n (t) — v Vu*(t) weakly in i 2 (ri;R JV ), we deduce that 

(»fe» + P?„ (*))(! - XsJ v «e„(0 ^ V«*(t) weakly in L 2 (ft;K w ) . (6.12) 

On the other hand, from the a priori estimate (|6.6[) . the Cauchy-Schwarz inequality, and (|6.9p . we infer that for 
every $ g %f c °°(ft; K^), 

1/2 / „ \ 1/2 



< 



< 



II*IIl~(0; R «) ^(1 + ?l £n )J? N (E n ) -> 0. (6.13) 



By (|6.6p and the boundedness of p £n , the sequence {(?? e „ + p £rl (t) 2 )\7u £n (t)} is thus bounded in L 2 (il;M. N ), so 
that (|6T2|) and (j6J3| yield 

(r? £n + p 2 (i))Vu e „(t) Vu,(*) weakly in L 2 ^;!^) . (6.14) 
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Finally, Lemmas 14.91 and 14.111 ensure that, for each n G N, the mappings t h-> (ij £n + p £n (t) 2 )Vu £n (t) arc 
strongly measurable in L 2 (il;M. N ). Hence t M> Vu*(i) is weakly measurable in L 2 (f2;]R Ar ), and thus strongly 
measurable owing to Pettis Theorem. The fact that Vu, G L°°(0, +oo; L 2 (ft; R N )) is a consequence of the 
second relation in (|6.3[) together with the uniform bound (|6.6[) . □ 



Our next goal is to pass to the limit as e n — > in the inhomogeneous heat equation solved by u £n . 

Proposition 6.3. The function solves the generalized heat equation 

'< = div(Vu*) - /3(it* - g) in L 2 (0, +oo; L 2 (Q.)), 
Vu«-^ = in L 2 (0,+oo;H- 1 / 2 (dCl)), 

u*(0) = u . 

Proof. By Corollary [IT?] u e n 1S the solution of the equation 

div((rfe M + P ljVu £n ) -/3(u en -g) in L 2 (0, +co; L 2 (fl)) 
= in L 2 (0, +00;^ / 2 (dfl)) 



du £n 
dv 



According to Proposition 16.21 u' £ + f3(u £n — g) — ^ u'^ + — g) weakly in L 2 OC ([0, +00); L 2 (S1)), so that it 
remains to pass to the limit in the divergence term. Let us fix $ G L 2 (0, +00; L 2 (£l;M. N )) arbitrary. Thanks to 
(|6.14[) . for a.e. t > we have 



(j? e „ + p 2 (t))Vu e „(t) -$(t)da; 



Vu*(i) • $(t)<fx. 



dt. 



By the dominated convergence theorem, we deduce that for every T > 0, 

/ / (Ve + pl )Vu e ■<&dxdt — ► / / Vu*-$dx 
Jo Jn " n ^°°Jo Jn 

Hence (?/ £n + p 2 n )Vu £n — ^ Vm* weakly in £ 2 OC ([0, +00); L 2 (f2)). On the other hand, it follows from the equation 
that the sequence {div((ry £ii + p £n )Vu £n )} is uniformly bounded in Lf oc ([0, +00); L 2 (il)). Up to a subsequence, 
it is therefore weakly convergent in L 2 OC ([0, +00); L 2 (f2)). Identifying the limits in the distributional sense, we 
infer that div((r/ £rl + p 2 n )\7u £n ) div (Vit*) weakly in L 2 QC ([0, +00); L 2 (Ct)) which yields 

< = div(Vw») - 0(u* - g) in L 2 {0, +00; L 2 (tt)). 



(6.15) 



Concerning the Neumann boundary condition, we first notice that Vu» G L°°(0, +00; L 2 (f2; R^)) and 
div(Vu*) £ L 2 OC ([0, +00); L 2 (f2)) imply that V«*(t) has a well defined weak normal trace for a.e. t > on <9f2 
(denoted by Vu*(i) • f), and that t h-> Vu*(i) • f belongs to L 2 OC ([0, +00); H~ 1 / 2 (dVL)). Then, given T > 0, we 
infer from the equation solved by u £n that 

/ / (u' Sn 0+ (rj E „ + p 2 jVu £ „ • V(j) + P{u £n -g)4>) dxdt = for all G L 2 (0, T; if^fl)) . 

Passing to the limit as e n — > in this variational equality yields 
i-t r 

«0 + Vu, • V0 + /3(w* -5)0) dxdt = for all G L 2 (0,T; H 1 ^)) 



Jo 



Together with l|6.15[) this last property ensures that Vu* • v = in i 2 (0, T; H~ 1 ^ 2 (dfl)) for every T > 0. 

Finally, the initial condition u* (0) = uq is a consequence of the fact that u £n (0) = «o together with the 
strong convergence in L 2 (il) of u £ „(0) to u*(0). □ 
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6.2. Limiting crack set and the energy inequality 

Our main goal is now to pass to the limit as e n — > in the energy inequality established in Corollary 15.151 We 
first notice that Theorem 12.21 and Proposition 16.21 (with A = 0) immediately imply that for every t > 0, 

= £ / \Vu*(t)\ 2 dx + ^ N - l {J Ur[t) ) + § / (u*{t)-g) 2 dx<]immfE Sn ( Uen (t),p Sn (t)). 

We emphasize that this lower bound only involves the measure of the jump set of u*(t). It will be later improved 
in Proposition 16 . 71 by replacing J u *{t) by a countably I ^f ,JV_1 -rcctifiable set T(t) containing J u ,{t) an d increasing 
with respect to t. This monotonicity property of the crack acts as a memory of the irreversibility of the process 
characterized by the non-increasing property of £ i — > p £n (t) together with the non-decreasing property of the 
diffuse surface energy © En established in Proposition 15.31 

To prove the assertion above, we fix an arbitrary countable dense subset D of [0,+oo), and we consider for 
each t G D and neN the bounded Radon measure 

Mn(f) := |Vp e „(i)r + - Pe n (t)) P ) J? N Ltl. 

V P P £n / 

By the energy inequality (|6.4j) . we infer that the sequences {p* n (t)} n eN are uniformly bounded with respect to 
t G D. Then, a standard diagonalization procedure together with the metrizability of bounded subsets of ^(R N ) 
yields the existence of a subsequence (not relabeled) and a family of bounded non negative Radon measures 
{/i(t)}t 6 D (supported in Q) such that 

fJ>n(t) n(t) weakly* in (R n ) for every t E D . 

We first claim that the mapping t G D u(t) inherits the increase monotonicity of the diffuse surface 
energy. 

Lemma 6.4. For every s and t G D with < s < t we have 

fi(s) < fl(t) . 

Proof. Let us fix s and t G D with < s < t. Let B C be an arbitrary Borel set, and K C B C A where 
A is open and K is compact. Let us consider a cut-off function £ G c tf£°(SL N ; [0, 1]) such that ( = 1 on K and 
C = on R N \ A, and let us define 

Pn :=0fe»(*) + (l-C)p e „(*)- 

Note that p n G W 1 '?^), and since t > s, we have p n < Pe„(s) in Q. As a consequence of the minimality 
property established in Proposition 15.31 we have 

/ (*£l\V Pen (8)\* + -?-(l-p en (8))A dx< [ (^-\Vp n \P +_2_(i_^ n )P > ) dx. 

JO V P P £ n ) Jn\ P p'Sn J 

Since Vp„ = (Vp En (t) + (1 — ()\7p En (s) + (p En (t) — p En (s)) VC, there exists a constant C > (independent of n) 
such that 

/ |Vpnpd!B< I |CVp e „(t) + (l-C)Vp e „( S )|P 

+ C / IVCKpeJ^-PeJ^ll + IV^Wr' + IVp^WI^ + IVCr^PeJO-PeJ^r 1 ) tfa 

Jn 

< jf (C|V Pen (t)|* + (1 - C)|Vp En ( S )r) dx + C(l + llVp^Wll^fn^j + ||Vp e „( S )||^ ;RjV) ) , 
where we used Holder's inequality and the fact that < p En < 1. Hence, 

Pn(s)(R N )< [ (dp n (t)+ [ (l-Odpn^ + CeP- 1 

+ Ceir 1]/ n\\^r 1)/p ^PeAmV {im + ll4 p - 1)/p Vp e ^)ll^ ;KN) ), (6.16) 

and passing to the limit as n — > oo yields 

p(s)(R N )<[ (du(t)- [ (du(s) + u(s)(R N ), 

JR N Jffi N 
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From this inequality we deduce that 



Taking the supremum among all compact sets K C B, the infimum among all open sets Ad B, and using the 
outer-inner regularity of the measures /j,(s) and fi(t) leads to /i(s)(B) < /i(t)(B). □ 

We can now define a family of increasing cracks for times in the countable dense set D. 
Lemma 6.5. There exists a family of countably J$? N ~ 1 -rectifiable subsets {T(t)} t £D of such that 

• r(s) C r(i) for every s < t with s, t G D; 

• J u (s) C r(i) for every t G D and < s < t; 



fi(t) > .jV N - l L-T{t) for every t G D. 



Proof. For t G D, let us define the upper density of /i(t) at x by 

O it, x) := hmsup rr— r- = lim sup rr— ;- 

r-vO ^N-\r" 1 r-S-0 WAT-ir JV 1 

for all x G R w , and the Borel set 

K(t) := {xeM. N : @*{t,x) > l} C H. 



Note that the monotonicity property established in Lemma [6741 ensures that D 3 t h-> 0*(t, x) is non-decreasing 
for every x G O. Consequently, 

if (s) C for every < s < i with s,t £ D . 
Moreover, from standard properties of densities (see [3j Theorem 2.56]) we infer that for every t G £>, 

j^JV-i [_!<-(£) < M (*). (6.17) 

Let us now fix i G D and s G [0,i] (not necessarily in D), and let A and A' C M be open sets such that 
A C A'. We consider a cut-off function C G ^(i"; [0, 1]) such that C = 1 on A n ft and ( = on R w \ A'. 
Arguing exactly as in the proof of Lemma 16.41 we obtain inequality (|6.16p from which we deduce that 

Hn(s)(A)< [ (d[i n (s)< [ (dvn(t) + Celr 1)/p , (6.18) 

for some constant C > independent of n. By (|6.5j) and (|6.3j) . we infer that 

\hnMii n (a)(A)>JF N - 1 (J u , a) riA). 

n—too 

Passing to the limit in f|6.18D then leads to 



Taking the infimum with respect to all open sets A' containing A yields 

ju(t)(I) > J^ ,N ~ 1 (J u ,(s) n A) for every open set A . 

In particular, since J u ,(s) is countably J^' JV_1 -rectifiable, we infer from the Besicovitch-Mastrand-Mattila The- 
orem (see [31 Theorem 2.63]) that Q*(t,x) > 1 for J0 s,JV_1 -a.e. x G J u ,(s), an d hence 

J„, (s ) C K(t) for every s G [0,t] . (6.19) 

The Borel sets {K{t)} t ^D have all the required properties, except that they might not be countably J*f N ~ 1 - 
rectifiable. However, since M' N ~ 1 (K{^)) < +oo by (|6.17|) . it is possible to decompose each K (t) into the union 
of a countably t ^ e,Ar_1 -rectifiable set T(t), and a purely J^' Ar_1 -unrectifiable set K(t)\T(t) (see e.g. [3J page 83]). 
This decomposition being unique up to Jff JV_1 -negligible sets, and J u ,( s ) being countably J^ Ar_1 -rectifiable, 
we deduce from (|6.19[) that 

J u ,(s) C f (t) for every t G D and s G [0, t] . 
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Moreover, for s,t G D with s < t we have T(s) C K(s) C K(t), and since T(s) is countably J^ 7V_1 -rectifiable 
we finally conclude that T(s) C T(t). □ 

We now extend our definition of crack set for arbitrary times. We set for each t > 0, 

r(t):= f| f(r). 

t>(,t£D 

Lemma 6.6. For every t > 0, the set T(t) is countably Jf 1 ^^ 1 -rectifiable, and it satisfies 

• r(s) C r(f) /or every < s < t; 

• J u , (t ) C T(t) for every i > 0. 

Proof. Clearly {r(t)} t >o is a family of countably Jff Ar ~ 1 -rectifiable sets satisfying T(s) C T(t) for every 
< s < t. Moreover, for t > we have 

Jf N - 1 (J Mt) \T(t)) = J? N -i (j u . (t) \ P| f(r) 

y T>t,reD 

= ( (J ( J «.(*) \ r(r)) I < £ (j Mt) \ f(r)) = , 

\T>t,TeD J T>t,T£D 

since J u ,(t) C r(r) for all t G £> such that r > t by Lemma |6"31 Consequently, J„,(t) C T(t). □ 

We are now in position to improve the energy inequality by replacing the jump set of u*(t) by the increasing 
family of cracks T(t) constructed before. 

Proposition 6.7. For every t > 0, 

\Vu4t)\"dx + .Jf N ~ 1 (T(t)) + ^ j {uS)-g?dx+ f \\<{s)\\l 2{n) ds 



<\f |Vu | 2 cfe + ! f (u -g) 2 dx. 
Proof. /Step 1. We first consider the case t £ D. According to the energy inequality in Corollary 15.151 we have 
(Ve n + p 2 Sn (t))\Vu en (t)\ 2 dx + n n (t)(m N ) + ^ j ( Ue „(t) - g) 2 dx + f H n (s)\\l HQ) ds < £ £n (u ,p £ n ) . 



* Jn z Jn Jo 

Since p, n (t) weakly* in J'fl") and > JF N - 1 \_f(t) by Lemma we have 

lim inf p n (t) (W N ) > > ^^(f (*)) . 

n—>oc 

On the other hand the second inequality in (|6.3j) with A = yields 

liminf / (r) e + p, (t) 2 )\Vu e {t)\ 2 dx > / |Vu*(t)| 2 dx . 
Jn " Jji 



Therefore. 



\ [ \Vu*(t)\ 2 dx + jr N -\t(t)) + ?- [ (u*(t)-g)*dx + tlK^Wl^ ds 
1 Jn A Jn Jo 

< liminf {lj(rk n + p 8n (t) 2 )|Vu 8n (f)| 2 dx + p n (t)(R N ) + | J^(u £n (t) - gf dx + J K»||i S(n) ds\ 



Then, by the minimality property (|3.1[) of p^" , we have 



S £n (u ,p^)<S £n (u ,l) = ^^- f \Vu \ 2 dx + ^- f (u - g) 2 dx -+ \ f \\7u \ 2 dx + £ / ( 

* Jo 1 Jn 1 Jn * Jn 



("o - gf dx . 
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which leads to 

f \Vu*(t)\ 2 dx + .Jf N ~ 1 (f{t)) + ^ I {u*(t) - g) 2 dx + I \K(s)\ 2 L2{n) d S 



2 In' w ' v v " 2 



< 



i f \Vu \ 2 dx + ^ f (u Q -g) 2 dx. (6.20) 



S'iep We now extend the inequality above to the case where t > is arbitrary. In that case, there exists a 
sequence {tj} C D such that tj — > t with tj > t. By (|6 . 20[) we have 



sup {||u*(tj)|U- ( n) + ||Vu»(^)]|^(n !R ^) + ^f^" 1 (-4, (**))} < 



oo . 



.ye: 



since J u ,{t-) C f (tj) by Lemma 15751 On the other hand, 6 ^4C 2 (0, +oo; L 2 (0)) by Proposition 16. 2[ and thus 
u*(tj) — > «*(i) strongly in L 2 (S1). Applying Ambrosio's compactness Theorem (Theorems 4.7 and 4.8 in [3]), 
we deduce that Vu{tj) — Vu(t) weakly in L 2 (fi;R Ar ). Since r(t) C r(tj) for all j £ N, we finally conclude that 



i JjVu*(t)\ 2 dx + ,^ N -\T(t)) + | y («*(*) - ff) 2 + y ||«*Wlli»(n)' 

< Uminf ji jf |V U *fe)| 2 <fa + ^- x (rfe)) + § jf («.&) ~ S) 2 + jf ' K(*)||i» (n) ds 
which, in view of (|6.20[) , completes the proof of the energy inequality. □ 



6.3. Relation with the unilateral slope of the Mumford-Shah functional 

In [TH] the authors have introduced a new notion of unilateral slope for the Mumford-Shah functional (in 
dimension N = 2 and for (3 = 0). In this paper, the Mumford-Shah functional is slightly different from the one 
we consider here (see (|2.3[) ) . It is rather defined on pairs (u, if) by 

£,(u,K) :=i f \\Ju\ 2 dx + .^ N - 1 {K) + ^- f (u - g) 2 dx , 
2 in 2 Jn 

where u £ SBV 2 (£l) and if is a subset of 51 satisfying J^ pN ~ 1 (K) < oo and J u C if. The related unilateral 
slope of £ * is then given by 

\o£* | (u, A ) := hm sup r n , 

«->u ||« - w|U=(n) 

where v — > u in i 2 (51). In [19], the authors proved that if |d£*|(u,if) < oo, then div(Vu) <E L 2 (r2), and that a 
weak form of 

£ = onif 

holds, where ^ denotes a unit normal vector on if. They also obtained the inequality |<9£*|(w, if) > ||div(Vu) — 
(3{u — g)\\L 2 (n), and that equality holds if u and if are smooth enough. By means of an explicit counterexample, 
they have shown that \d£J\ is not lower scmicontinuous for any reasonable notion of convergence. In view of 
this result, they have introduced a notion of relaxed slope corresponding to a lower scmicontinuous envelope of 
<9£» with respect to a suitable sequential topology. More precisely, the relaxed slope \d£* \ is defined for a pair 
(u, if) in the domain of by 

|<9£* | (it, if) := inf { lim inf \d£ * | (u„ , if „) } 

where the infimum is taken over all sequences {(u n , if n )} n eN such that u n — > u strongly in L 2 (S1), Vu„ — ^ V« 
weakly in L 2 (fl; R N ), and if„ cr 2 -converges to if (see [14] Definition 4.1] for a precise definition). They established 
that if |9£»|(ti, if) < oo, then there exists / 6 L 2 (fl) such that 

' -div(Vu) = / in L 2 (fl), 

|Vm| 2 - div(itVu) < fu m^'(fi), (6.21) 
Vii'i/ = inii- 1 / 2 ^). 
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Again, there is an inequality \d£*\(u, K) > ||div(Vit) — j3{u — g)||z, 2 (fi)j and equality holds in some particular 
cases. Note that, in the case where u and K are smooth enough, the first line in (|6.2f j) implies the continuity 
of jf£ across K, and the second one is then a weak reformulation of 

du 

(u + -u~)— >0 onK, 
ov 

where are the one-sided traces of u on K according to the orientation v. 

In our context, the analogy between the definitions of the unilateral slopes \d£ E \ and \d£*\ is quite clear, and 
it was actually one of the motivations to introduce \d£ E \. In view of the relation between the Ambrosio-Tortorelli 
functional and the Mumford-Shah functional in terms of T-convergence, a very interesting issue would be to 
find a precise relation between \d£*\ and the asymptotic behavior as e \. of \d£ E \. Even if we do not pursue 
this issue here, we prove for completeness that a conclusion similar to [TSJ Proposition 1.3] holds for \d£ E \. For 
clarity reasons we only state the result in terms of the asymptotic limit u* previously obtained. 

Proposition 6.8. For a.e. t> 0, we have 

||div(Vu.(t))-0(u,(f)-$)|| La (n) < liminf \d£ En \(u En (t), p En (t)) < oo , (6.22) 

v ' n— too 

and 

|Vw*(t)| 2 -div(u*(i)Vit*(i)) < div(Vu,(t)) in 2>'{tt). 

Proof. From Propositions 13.81 fc 15.141 together with (|6.4[) and Fatou's lemma, we first deduce that 

r + °o r+ca 

/ liminf \d£ e \ 2 {u £ (t),p £ (t)) dt < liminf / \d£ e \ 2 {u £ (t),p e (t))dt<C, 
Jo n ^>°° ™^°° Jo 

for a constant C > independent of n. Hence there exists an Jz? 1 -negligible set C C (0, +00) such that 

liminf \d£ En \ (u en (t), p En (t)) < 00 for t £ (0, +00) \ C . 

n— >oc 

Let us now fix t S (0, +00) \ C and extract a subsequence (depending on t) such that 

lim \d£ en .\(u en .(t),p Sn .(t)) = liminf \d£ erl \(u 6n (t), p £n (t)) . 

j— >00 3 3 3 n — j-oo 

By Proposition 13.91 the sequence {div((?7 e? ^ + p\ n (t))Vu En . (t)) } is thus bounded in L 2 (Q), and in view of 
(I6.14p we deduce that 

div((r/ £nj + p 2 en {t))Vu £nj (t)) div(Vu*(i)) weakly in L 2 (Q,) . (6.23) 

Then (|6.22|) follows from the convergences in (|6.1|) . and the lower semicontinuity of the L 2 (J7)-norm. 
Using again Proposition 13.91 we next notice that 



(Ve n + p 2 eJt))\Vu £n (t)\ 2 <pdx + j U en (t)(n £n + p 2 En (t))Vu En (t) ■ Vipdx 

u En {t) div((?7 e „ + p 2 £ri {i))Vu En {t)) ipdx 
for any nonnegative function ip € 2$(Q), and the conclusion follows from (|6.14[) and (|6.23[) . □ 
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